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ABSTRACT 

A unified description is provided for the 21cm signatures arising from minihaloes 
against a bright background radio source and against the Cosmic Microwave Back- 
ground (CMB), within the context of a dynamical collapsing cosmological spherical 
halo. The effects of gas cooling via radiative atomic and molecular processes and of 
. star formation on setting the maximum mass of the minihaloes giving rise to a 21cm 

• ' signal are included. Models are computed both with and without molecular hydrogen 

\ formation, allowing for its possible suppression by an ambient ultra-violet radiation 

field. The spectral signatures and equivalent width distributions are computed for a 
Q ' ACDM cosmology. The detectability of minihaloes in absorption against bright back- 

, ground radio sources is discussed in the context of future measurements by a Square 

c/3 ■ Kilometre Array (SKA) and the LOw Frequency ARray (LOFAR). The brightness 

I temperature differential relative to the CMB is also computed. 

Several generic scenarios are considered. For the cosmological parameter con- 
' straints from the Wilkinson Microwave Anisotropy Probe {WMAP)^ in the absence 

^ , of any form of galactic feedback, the number of systems per unit redshift in absorp- 

psj ■ tion against a bright radio source at 8 < z < 10 is dN/dz ~ 10 for observed equivalent 

^sO ' widths exceeding 0.1 kHz. For larger equivalent widths, somewhat fewer systems are 

CO , predicted at increasing redshifts. The estimated numbers are independent of the pres- 

ence of star formation in the haloes following molecular hydrogen formation except 
I for rare, high equivalent width systems, which become fewer. LOFAR could plausi- 

, bly detect a minihalo signal against a 30 mJy source in a 1200 hr integration. SKA 

could detect the signal against a weaker 6 mJy source in as little as 24 hrs. Adding 
cosmological constraints from the Atacama Cosmology Telescope (ACT) suppresses 
the predicted number of all the absorbers by as much as an order of magnitude. In 
the presence of a background of ambient Lya photons of sufficient intensity to couple 
the gas spin temperature to the kinetic temperature, as may be produced by the first 
5^ ■ star-forming objects, the number of weak absorption systems is substantially boosted, 

by more than two orders of magnitude, rendering the signal readily detectable. Weak 
absorption features arising from the cold infalling regions around the minihaloes may 
appear as mock emission lines relative to the suppressed continuum level. A moderate 
amount of heating of the IGM, however, would greatly reduce the overall number of 
absorption systems. 

By contrast, the absorption signal of minihaloes against the CMB is distinguish- 
able from the diffuse IGM signature only for a limited scenario of essentially no feed- 
back and moderate redshifts, z < 19. The strength of the signal is dominated by the 
more massive minihaloes, and so sensitive to a cutoff in the upper minihalo mass range 
imposed by any star formation and its consequences. Once the first star-forming sys- 
tems provide feedback in the form of Lya photons, the diffuse IGM signal will quickly 
dominate the signal from minihaloes because of the small total fraction of IGM mass 
in the minihalo cores. 

Key words: cosmology: theory - dark ages, reionization, first stars - galaxies: for- 
mation - intergalactic medium - molecular processes - radio lines: general 
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1 INTRODUCTION 

A principal science goal of the Square Kilometre Array 
(SKAj3 is the detection of neutral hydrogen in the Inter- 
galactic Medium (IGM) following the recombination epoch 
and prior to the Epoch of Reionization (EoR), performing 
3D tomography of the filamentary large scale structure of 
the Universe as revealed by neutral hydrogen at the end 
of the Dark Ages when the first so urces of light emerge 
in the Universe after the Big Bang (|Hogan fc ReesI 1 1979 : 
IScott fc Re3ll990l : lMadau et al.lll997l ). Even before SKA is 
built, the absorption of a still largely neutral IGM may be 
detectable against radio loud galaxies or quasars by SKA 
pathfinders like the LOw Frequency ARray (LOFARjfl and 
the Long Wavelength Array (LWAjf], both already under- 
way, and the Murchison Widefield Array (MWAfl, a SKA 
precursor facility soon to start. 

Extracting the EoR signal from the radio data is a 
major undertaking requiring modelling of both the ex- 
pected 21cm signature and the effects of contaminating 
sources which will swam p the underl ying cosmolog ical sig- 
nal LHogan fc Roes 19Zi; IShaver et al. 1999; Moral es etal] 
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l2006l : ijelic et al.|[200^ '). The extraction of the EoR signature 
is further complicated by uncertainties in the signature it- 
self arising from the unknown manner in which the Universe 
was reionized and from the structure of the IGM. The 21cm 
signals will be highly patchy, as the cold hydrogen atoms 
will closely trace the filamentary structure of the dark mat- 
ter in the Universe on large scales and the non-linear fluc- 
tuations in the dark matter on the very small. If a bright 
background radio source is present, the intervening neutral 
hydrogen in the IGM will produce an absorption signature 
against it, with a strength depending on the temperature of 
the I GM and how efficiently the spin temperature is coupled 
to it l|Fieldlll959"al '). The IGM will also produce a fluctuating 
signal ag ainst the Cosmic Microwave Background (CMB) 
l|Hogan fc Ree s 1979; Scott & Rccs,1990; Tozzi et al. 200a), 
either in absorption or emission depending on the gas kinetic 
temperature compared with the CMB temperature. Only on 
scales small compared with the Jeans length of the gas will 
the IGM produce a smooth signal. The neutral hydrogen 
thus provides a unique means of measuring the power spec- 
trum of dark matter over a wide range of length scales span- 
ning six decades at the as yet unprobed epochs between the 
Recombination Era and the emergence of the first galaxies. 

Perhaps the best near-term prospect for detecting a still 
neutral IGM is in absorption against a bright background 
radio source. Collapsed minihaloes will produce a 21-cm for- 
est of absorption features in the radio spectrum of such a 
source (|Carilli et a l. 2002; Furlanctto & Locb 2002). Whilst 
this will be superposed on a global absorption signature aris- 
ing from the fluctuating diffuse IGM, the deeper absorption 
features may be more amenable to a clean detection than the 
signal against the CMB by an instrument able to resolve the 
minihalo signature in frequency. 

Two types of experiments may detect the 21cm signa- 
ture against the CMB, a direct brightness temperature dif- 
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ferential based on differencing the signals between a neu- 
tral patch and an ionized patch, and fluctuations in the 
brightness temperature based on differencing the signals be- 
tween two (or more) neutral patches separated in either fre- 
quency or angle. The detection of either signature requires 
that the hyperfine spin structure of the hydrogen be decou- 
pled from the CMB. The signal will consequently reflect the 
processes which achieve the decoupling in addition to the 
spatial distribution of the gas. Two mechanisms dominate: 
decoupling through collisions with electrons and other hy- 
drogen atoms, and decoupling resulting from the scattering 
of Lyof (and other Lyman re sonance line) photons through 
the Wouthuysen-Field effect (|Wouthuvsenll 19521 : lFieldlll958l ) 
once a sufficient intensity of Lyman resonance line photons 
builds up in the IGM from the first stars and galaxies. 

Even without an adequate supply of Lyman resonance 
line photons, collisions provide a sufficient means of cou- 
pling the spin temperature to the gas tempera ture in struc- 
tures with overdensities 5p/p > 30 a t a > 6 l|Madau et al.l 
1 19971) . It was argued by llliev et all (|2002D on the basis of 
semi-analytic modelling that minihaloes will produce a sub- 
stantial signal against the CMB at z > 6, and will dominate 
over the diffuse IGM before an adequate flux of Lya pho- 
tons develops to decouple the spin temperature of the dif- 
fuse IGM from the CMB temperature. A differential bright- 
ness temperature based on differencing between two neutral 
patches will in particular be enhanced b y halo bias, as well 
as th e bias of the reionization sources (iMorales fc Wvithd 
I2OIOI ). This claim has been disputed bv lOh fc Mackl (|2003l ). 
who objected that, because of the overwhelming amount of 
mass in the diffuse IGM compared with minihaloes, except 
very early on before there is any feedback from galaxies, 
the diffuse component will dominate. It has been argued 
bias boosts the brightness temperature fluctuations signif- 
cantly only on scales too sma ll to be measured with cur- 
rently planned radio facilities (|Furlanetto fc Ohll2006l ). Nu- 
merical simulations have tended t o support these expecta- 
tions (Kuhlcn ct al. 2006i: lYue et al. 2009). Usi ng an exten- 
sion o f the Press-Schechter formalism, iFurlanetto fc Loebl 
(j2004 ) suggested that even without galactic feedback, the 
21cm signal from the shocked diffuse IGM would exceed 
that due to the minihaloes. H igh spatial resolution numerical 
simulations bv iKuhlen et al.l (2006) appear to confirm this 
cla im, although the opp osing view has again been defended 
bv lShapiro et al.l l|20od ). 

Estimates of the strength of the 21cm signature from 
minihaloes are made uncertain by several poorly known fac- 
tors. The number density of the haloes is notoriously difficult 
to quantify, especially since the minihaloes arise on scales for 
which the dark matter power spectrum is approaching the 
"flicker noise" small-scale limit, for which structures on all 
scales collapse and virialise simultaneously. Numerical sim- 
ulations are currently not of much help for estimating the 
distribution of dark matter haloes or their internal structure 
on all the relevant scales, as it is not yet possible to capture 
the power over the full dynamic range required. 

Another major outstanding uncertainty is the effect of 
the first generation of stars on subsequent star formation. 
Early models suggested Lyman- Werner photons from 
the first stars will dissociate molecular hydrogen in the 
surrounding gas and inhibit star formation in haloes with 
temperatures below 10'' K, when atomic cooling dominates 
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l|Haiman et all Il997al lbh. Subsequent analyses have shown 
the picture is much more complicated. Self-shielding will 
limit the radiative transport of the Lyman- Werner pho- 
tons, although by an amount that ge nerally depends on 
internal flows (jClover fc Brandl 1200 ll ). Whilst hot stars 
will photoionise their surroundings, partial ionisation may 
enhance the rate of molecular hydrogen formation, as 
would sources of x-rays, negating the disso ciating effects 
of th e UV radiation from early sources l|Haiman et al.l 
|2000| ). Within relic H ii regions, molecular hydrogen 
reforms as the gas recombines, re-enabling star formation. 
The injection of dust would provide still an additional 
mechanism for molecular hydrogen formation. The injection 
of metals from supernovae would accelerate cooling. The 
feedback from star formation on subsequent generations 
of st ars has been cons i dered widely in the lite r ature , 
eg. jCouchnian fc Reej 1 19861: iMac Low fc Shulj 1 19861: 



1997al: 



Dekel fc Reeslll987l: lOstriker fc Gnedinlll996l: iHaiman et aD 



Omukai fc Nishi '199^; 'Glover & Brand' '2001'; 
Machacek et al., 2001: O'Shea ct al. 2005; Yoshida et al., 
20071: IWise fc Abel l2007l: Iwhalen et ahl I2OOS 



Mesinger et al.ll2009l : IWhalen et al.H2010l : lGnedinlf2010l ~ 

Even if the first stars produce local pockets of Lyman- 
Werner photons, it is unknown whether or not a metagalac- 
tic radiation field ever develops that is sufficiently intense to 
dissociate molecular hydrogen in minihaloes everywhere, or 
the epoch at which such a field develops if it does. Indeed, 
one of the primary goals of the search for an intergalactic 
21cm signature is to establish the epoch of star formation 
and quantify its large scale distribution. To allow for the 
possibility that the first generation of stars either only par- 
tially or completely inhibits subsequent molecular hydrogen 
formation, two sets of models are computed, one allowing for 
molecular hydrogen formation and a second with molecular 
hydrogen formation suppressed. 

Early massive stars may also produce x-rays, par- 
ticularly after they explode as supernovae and heat 
their surroundings to high temperatures. Other possible 
sources of x-rays include early Active Galactic Nuclei and 
collapsed structures shock-heated to x-ray temperatures 
jMadau et al.l Il997): iTozzi et al.l I2OO0I: ICarilli et al.l |2002| : 



iFurlanetto fc Loeb|]2004l : iKuhlen et a l. '2006"). The amount 
of heating is unknown. Substantial heating would reduce the 
collapse fraction of the gas into minihaloes, weakening their 
21cm signature. But even a small amount too little to much 
affect the dynamics of the gas will reduce the absorption 
signal against a bright background radio source, and could 
change an absorption signal against the CMB into emission 
if the IGM temperature increases above that of the CMB. 

The primary purpose of this paper is to address the 
effect of some of these assumptions and uncertainties on 
the expected 21cm signatures from minihaloes. Only single- 
point statistics are addressed, relating to the absorption 
signature of minihaloes against a bright background radio 
source, or the differential brightness temperature against 
the CMB based on comparing the signal along the line of 
sight through a neutral patch with that through an ionized 
patch. Cross-correlation measurements, particularly as have 
been considered in the context of measurements against the 
CMB, require estimates of the halo bias factors and spatial 
correlations between the minihaloes. Correlated fluctuations 
arising from the diffuse IGM are also expected to be a ma- 



jor contributor. These topics are beyond the scope of this 
paper. 

In previous models of the 21cm signature from mini- 
haloes, the minihaloes have generally been assumed to be 
isothermal and in hydrostatic equilibrium. It was recognised 
in the context of the minihalo model for Lya forest ab- 
sorption systems ijlkeuchil Il986l : iReed 1198^ ). that the sys- 
tems will iri gene ral be formed in dy namically colla psing 
iBond et al.lll988l ') and non-isothermal (iMeiksinlliggj ) dark 
matter haloes. Allowing for the dynamics introduces mod- 
ifications to the thermal structure and pressure support of 
the gas in the minihaloes, which in turn modify the cross 
sections for their detection. 

The haloes are evolved as spherical density perturba- 
tions. Treating the haloes as arising from isolated perturba- 
tions is itself an approximation, as in bottom-up scenarios 
like cold dark matter dominated cosmologies, structures are 
created from the merger of smaller systems. In this sense, the 
perturbations are more representative of an evolving patch 
of smaller haloes which merge into the final collapsed struc- 
ture. In the early stages of evolution, even this is not a good 
approximation as the overdense patches are typically triaxial 
(|Bardeen et al.|[l986l ) (hereafter BBKS), and em bedded in a 
large-scale filamentary network, the cosmic web ("Bo nd et al.l 
I199Q1 . Simulations suggest, however, that spheroidal struc- 
tures do form; the spherical minihalo model provides a useful 
description of the higher column density systems within the 
Lya forest (|Rauchlll99i lMeiksinll2"009h . 

The collapsing spherical halo approximation has the ad- 
vantage over the static isothermal halo approximation of 
allowing the density and temperature profiles to arise nat- 
urally during the collapse of the halo. The modelling of an 
evolving system also reveals some of the dynamical effects 
which may play an important role in the cooling of the gas 
and in determining the minimum halo mass at which star 
formation sets in. It is argued that only a small number of 
stars need form before the remaining gas in the minihalo 
will be expelled through mechanical feedback due to photo- 
evaporative winds and supernovae. The rate of star forma- 
tion is modelled to provide an estimate for the upper limiting 
minihalo mass at any given collapse epoch. (More massive 
systems than minihaloes become increasingly complex, such 
as forming galaxy haloes and their associated discs, and are 
not considered here.) 

Only the 21cm signal from fully collapsed haloes, along 
with their surrounding infall regions, is considered in this 
paper. Sufficiently dense gas within the filaments will also 
contribute to the 21cm signal. The treatment of such mod- 
erate overdense structures is less straightforward to model 
and so is deferred to a later work, although comparison is 
made with the total signal that could potentially arise from 
a diffuse homogeneous component. Ultimately a very high 
resolution fully 3D combined dark matter and hydrodynam- 
ics simulation would be required for a precise prediction of 
the expected 21cm signatures. 

A second purpose of this paper is to provide a semi- 
analytic framework that incorporates the dynamical effects 
to estimate the 21cm signals from minihaloes. Numerical 
simulations are still unable to resolve the full range of scales 
required to estimate the 21cm signature without some semi- 
analytic modelling, and even if they were able to do so, they 
are very costly to run. Alternative halo mass distributions 
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may be readily explored within the semi-analytic framework. 
The 21cm signatures produced by minihaloes probe a part of 
the cosmological dark matter power spectrum that has never 
been measured at high redshifts, down to scales on the order 
of a comoving kiloparsec, or a millionth of the cosmological 
horizon. The effect of different cosmological scenarios, such 
as alterations to the small-scale power spectrum, may be 
readily examined within the framework. It is demonstrated 
below that the minihalo 21cm signature may in fact provide 
a useful means for distinguishing between rival predictions 
for the amount of power on small scales. 

Unless stated otherwise, the present-day values of the 
cosmological parameters assumed are flu = 0.27, Qv = 
0.73, Qb = 0.0456, h = 0.70, crg^-i = 0.81 and n = 0.96 
for the total mass, vacuum energy and baryon density pa- 
rameters, the Hubble constant (ft = _ffo/100kms~^ Mpc~^), 
the linear density fluctuation amplitude on a scale of 
8h~^ Mpc and the spectral index, respectively, consistent 
with CMB measurements by the Wilki nson Microwave 
Anisotropy Probe (WMAP) (|Komatsu et al. 2009. 2011). 

This paper is organised as follows. The next section 
presents basic relations for evaluating the 21cm signatures 
from minihaloes, and applies these to the simple model of 
tophat minihaloes with uniform densities that provides a 
fiducial for comparing the more complete dynamical models 
against. Sec. [3] describes the dynamical minihalo model and 
presents the 21cm signatures arising from individual haloes. 
In Sec. |4l the individual minihalo signatures are combined 
to model the expected statistical signals from the haloes. 
The statistical signals are discussed in various cosmological 
scenarios in Sec. (5] A summary of the principal findings and 
conclusions is provided in Sec. [S] Semi-analytic estimates 
for the post-collapse temperatures of minihaloes, their abun- 
dances, and the computation of molecular hydrogen forma- 
tion are discussed in separate appendices. 



2 21CM SIGNATURE OF TOPHAT 
MINIHALOES 

2.1 Spectral characteristics of top- hat minihaloes 

The diffuse component of the IGM at redshift z will pro- 
duce a relative differential brightness temperature against a 
background source of observed antenna temperature Tb{0) 
of 



STb{0) 
Tb{0) 



Ts 



Ts(0)(l 
1 - 



[1 - exp(-rd) 



Ts 



Tb{0){1 + z) 



(1) 



where the 21cm optical depth of the IGM with mean hydro- 
gen density nH(z), neutral fraction xm{z), and spin temper- 
ature Ts{z) is given by 

= J^^Hifl + 5{z)]nu{z)Xlo yj*^^ 



0.0046[1 + S{z) 
1 



( xiii{z 



[dv{l)/dl]/H{z) 
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-1/2 



{1 + Z 



,3/2 



1 - D,„ 



1 -1/2 



nm{i + z) 
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with Tb = rs(0)(l -I- z), and where T* = huio/ks, 
Aio — 21.1 cm, vio — c/Xio, fee is the Boltzmann constant, 
Aio — 2.85 X 10~^^ s^^ is the spontaneous transition rate of 
the 21cm hyperfine transition, and allowing for linear den- 
sity fluctuati ons S(z) and total line-of-sight velocity gradi- 
ent d v(l)/dl (Field 'l959al: IScott fc Reej|l990l : iMadau et all 
ll997l : iMeiksin.2009 ). The spin temperature is given by 



Ts 



T'cmb(2) + yaTa + VcTk 



(3) 



1 + ya + yc 

("Field! Il958l ). where y„ = (PioMio) (T./T^) for a de- 
excitation rate Pio of the hyperfine triplet state by Lya pho- 
ton scattering at the rate Pa ~ (27/4)Pio, ya is the colour 
temperature of the radiation field, yc — (Cio/j4io)(T«/Tk) 
for a coUisional de-excitation rate Cio ~ ki-oWhi of the 
hyperfine triplet state, where ki_o is the de-excitation 
rate coefficient, Tk is the kinetic temperature of the gas, 
and Tcmb{z) is the CMB temperature at redshift z. The 
colour temperature relaxes rapidly to the gas kinetic tem- 
perature after 10'' scattering times (Field 1959b; Mciksi^ 
I2006I ). so that Ta = Tk is assumed. For a Lya scat- 
tering rate Pa matching the thermalization rate Pth = 
(27/4)AioTcmb(z)/T», the radiation field effectively couples 
the spin temperature to the kinetic temperature of the gas 
( Madau et al. 1997,). The coUisional de-excitatio n rate co- 
efficie nt Ki_o is interpolate d from the tables of IZygelmaj] 
(|2005l ) for r < 300 K and lAllison fc Dalgarnd (| 19691 ) for 
T > 300 K, using 4 / 3 the latter's tabulated values as advo- 
cated bv lZveelmaiil (|2005l ). The optical depth will typically 
be much smaller than unity. The signal will be either in emis- 
sion or absorption against the background source, depending 
on the larger of Ts and Ts(0)(l -I- z). 

Superposed on any signature from the diffuse IGM in 
the spectrum of a background radio source will be narrow 
features produced by intervening minihaloes. These again 
will be either in emission or absorption, depending on the 
brightness temperature of the source. Whilst the collapsed 
haloes will develop a steep internal density profile, it is useful 
to formulate the contribution of the minihaloes in the sim- 
plified scenario of collapsing isothermal spherical tophats. 
Although a simplification, it provides a straightforward fidu- 
cial model that incorporates most of the essential physics 
against which more sophisticated models may be compared. 
In particular, in the limit of negligible dissipation, it pro- 
vides the expected scaling relations of the mean equivalent 
width on the cosmological parameters and the redshift of 
collapse. 

The differential brightness temperature produced by in- 
tervening gas against a background source of brightness tem- 
perature TBiy) is given by 



5Tb{v) = 



dlTs(l)'^e--"^'' 



TBiy) (1 -e" 



where 



— AioAoi 



d/'i2;Hl(i')"H(?')V'^ 



(4) 



(5) 



is the optical depth from the source up to a distance I 
through the intervening gas of total hydrogen density n-niV), 
neutral fraction xy{_i{V) and spin temperature Ts{l). The fac- 
tor 1/4 accounts for the occupation fraction of the lower 
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hyperfine level of the hydrogen. The factor T,/Ts accounts 
for stimulated emission of the 21cm line. The absorption 
line profile is described by ipi, = 4'v//S.ud, where <j)v = 
'K~^^'^ exp{ — [;^ — :/io(l — v(0/c)]^/(Az/d)^} is the dimension- 
less Doppler profile with Doppler width /S.ud = {h/c)u\o and 
Doppler parameter b — (2kBT/muy^'^ , where T is the gas 
temperature. Estimates for the gas temperature in collapsed 
minihaloes are provided in Appendix [X] The Doppler shift- 
ing of the line centre frequency by the bulk motion of the 
gas with fiow velocity v{l) along the line of sight has been 
allowed for in the optical depth. 

The observed equivalent width through an individual 
tophat halo at a projected comoving separation b± from the 
centre is given by 



obs / 7 



±(1 + 2) 



= ±il + z) 



dv 



STb{u) 
Tb 



Ts 
Tb{z) 



where the redshift factor converts the rest frame equiva- 
lent width to the observed frame, and the sign conven- 
tion is chosen so that the equivalent width is made posi- 
tive whether corresponding to emission or absorption. Here, 
= tq-k^^'^ Avd'Pv ■ Assuming the gas is isothermal and in 
hydrostatic equilibrium, the line centre optical depth tq for 
a spherical tophat collapsed at redshift Zc is given by 



To 



87r3/2 



. ,3 r-2Xmn'a{zc) 2ro 



4fe 



1 + Zc 



1 



fcTO 
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where the halo virialises at radius r„ = /cro/(l -I- Zc) with 
/c = {18n^)~^^^ . The factor accounts for the increased 
column density through the collapsed halo. 

Typically tq << 1, in which case the observed equiva- 
lent width of a halo of mass M collapsed at redshift Zc may 
be expressed as 
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fcTO 
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where 



<=(0) = 



(1 + z,m 
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^a;mnH(0)2ro(M) 
Ts(zc) ^ 
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Tb{z^) 

It will be shown below that for either absorption against a 
bright radio source or for absorption or emission against the 
CMB, the observed equivalent width is independent of the 
redshift of the minihalo, depending only on its mass. 

The combined absorption by the minihaloes will reduce 
the antenna temperature by exp(— r;), so that 5Tb{0) = 
-Tb(0)(1 - e"^'), where 

_{l + z,f 



±- 



1^10 



J obs ^ ^ obs 
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where N/dw^g^dz is the number of features of ob served 
equivalent width w"]^" per unit redshift l|Meiksinll20"09l '). Val- 
ues with either r; > or r; < are possible, corresponding to 
absorption or emission, respectively, against the radio back- 
ground. 

The equivalent width distribution may be computed as 

d'^N/dw^„dz = -{d/dw°^;)dN{> wZl'')/dz, where dN{> 



■w°^^)/dz is the number of absorbers per unit redshift with 
observed equivalent widths exceeding w°g^. This is given by 



dN{> <fa°) 
dz 



dip 
dz 



dlosM 



1 + z 



(11) 



where dlp/dz = {c/H{z)){l + z)~^ is the differential proper 
length per redshift, dn/dM is the number density of col- 
lapsed haloes of mass M per unit comoving volume, and 
o-th''(A'/) = 7r[6™'"'(M)]^ is the comoving cross-section cor- 
responding to the maximum comoving impact parameter 
bf'^^{M) through the halo within which the equivalent width 



exceeds From Eq. dSJ, 
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it is given by 

1/2 



(12) 



The last form expresses the number density per unit redshift 
more generally in terms of the comoving mean free path 



-^mfp(^ 



l/idn/dlogM)a^S^^{M), where 0\'i(M) is 



the cross section through a halo of mass M giving rise to an 
absorption feature with observed equivalent width exceeding 

It follows from Eq. pup that the cumulative optical 
depth from the ensemble of minihaloes may be expressed, 
after an integration by parts, as 
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where the last form is general, with E^J^^ ~ Uh{w[ 
equivalent-width weighted cross section, where {w°^ 
mean observed equivalent width face-averaged over a comov 
ing halo cross-section an- More generally. 



Eobs / I maxN 
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= 27r 



db±b±w2^^{bA 



(14) 



This may be compared with the optical depth through 
the diffuse component of the IGM. The spin temperature 
will quickly couple to the CMB temperature, in which case 
the optical depth of the diffuse component would Td — 0.011 
at z = 10. In the presence of either a sufficient collision 
rate with hydrogen atoms and electrons, or a sufficiently 
intense Lya scattering rate, the spin temperature of the 
hydrogen will be coupled to the kinetic temperature of the 
gas. The resulting optical depth of the diffuse component 
would then increase to ~ 0.13. In the presence of heating, 
the optical depth would be reduced inversely with increasing 
gas temperature. 

Experiments designed to measure an IGM signal against 
the CMB compare the signal from a neutral patch of the 
IGM with that towards an ionized patch, through which the 
CMB signal will be negligibly attenuated. For a spin temper- 
ature coupled to the gas kinetic temperature throughout the 
IGM, the corresponding differential brightness temperature 
compared with the CMB would be Tb — Tcmb — —330 mK 
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at z = 10. If the IGM were heated to a temperature much 
higher than the CMB temperature, the differential bright- 
ness temperature against the CMB, however, would convert 
to emission, saturating at a value of 29 mK. 

The minihalo signal against the CMB may be detectable 
when the amount of galactic feedback is too little to decou- 
ple the spin temperature of the diffuse IGM component from 
that of the CMB. Since the telescope beam probing the neu- 
tral patch will encompass a large number of minihaloes, the 
collective optical depth of an ensemble of minihaloes deter- 
mines the strength of the minihalo signal. The expression 
for the optical depth takes on a particularly simple form 
in the approximation that the systems are optically thin. If 
T(Af) denotes the hydrogen mass fraction of a halo of to- 
tal mass M, then using Eq. Q, the optical depth may be 
re-expressed as 

(1 + -)^ 



'32^^"^^« H{z) Tcmb(z) 



1 



dM-^T{M)Mr,cMB{M), 



(15) 



where 



VCMB 



Tcmb(zc) 



Ts 

rcMB(z) 



Tcmb{^) 



- 1 



l + #- 

-' th 



' rCMB(z) 



(16) 



is the 21cm efficiency for absorption (Tg < T'cmb(zc)) or 
emission (Ts > rcMB(zc)) against the CMB (|Madau et al.l 
Ll997). The last form uses Eq. ((3]). Comparison with 
Eq. ^ shows the optical depth takes on the same 
form as for the diffuse IGM with Ts{z) replaced by 
Tcmb(z) and nniz) by nn{z){fM), where (/a/) = 
[mHnH(0)]~^/ dM(dn/dM)r{M)MrjcMB{M). The optical 
depth is thus proportional to the mean hydrogen mass frac- 
tion of the IGM in haloes, weighted by the 21cm efficiency of 
the haloes. The resulting observed temperature differential 
compared with the CMB is 5Tb{0) ~ -r,rcMB(0). 



2.2 21cm signature against a bright radio source 

The diffuse IGM will produce an absorption signature 
against a bright background radio source such as a quasar 
or radio galaxy. For a mean neutral hydrogen density given 
by the cosmic mean hydrogen density, and for a hyper- 
fine structure coupled strongly to the CMB, so that Ts = 
TcMBiz), the optical depth of the diffuse component is 
Td — 0.011[(1 + z;)/ll]^^^. The absorption will appear as 
a step shortward of the observed wavelength of the 21cm 
line. 

Superposed on the step will be deeper absorption fea- 
tures arising from overdense structures, the 21-cm forest. 
The equivalent width of a minihalo with line centre optical 
To{b±) a projected distance b± from the cloud centre is 

,_i26 



(1 + ^) 



-//ioF[ro(6x)], 



(17) 



where F{to) is a function of the line centre optical depth. For 
a spin temperature coupled to the post-shock temperature 
of the halo, using Ea. HA4p from Appendix 1X1 in Eq.((7| gives 
for the line-centre optical depth through the core {b± = 0) 



ro(0) ~ 0.00275(1 + zc) 



1/2 



M 



-2/3 



(18) 



lO^M©^ 

Only the lowest mass haloes at high redshift will have 
ro(0) > 0.1, so that the linear curve-of-growth approxima- 
tion F(ro) ~ (7ri/V2)To may generally be made. 

The corresponding equivalent width may then be re- 
expressed as 

l2/3 / u^rj. X 1/3 



nH(0) 



S/cbT. 



0.2933 



0.0252 



47rGri™,pcrit(0)m 



M 



(- 



\GMfh ) 

1/3 

kHz 



M 



-1/3 



kHz. 



VlOSAfo; ^^^^ 

It follows that the observed equivalent width is independent 
of the collapse redshift Zc- The strongest absorption lines 
arise from the lowest mass haloes because of their lower 
temperatures. From the dependence of ro and w'^'^iQ) on 
M, and noting from Figure iBll that M^dn/dM — ^ constant 
for low mass haloes, it follows that the absorption signal 
varies like ti ~ A^mh/"^) so that it is the lowest mass haloes 
A/min, those just Sufficiently massive for shocks to form when 
they collapse, that dominate the overall signal. 

The characteristic observed line width of the features is 

1/3 

kHz. (20) 



Az.^'^^- 5.17(1 2c) 



-1/2 



/ M 
V W^Mq 



Such features would be well resolvable by an instrument with 
channels on the order of 1 kHz wide. 



2.3 21cm signature against the CMB 

In addition to an absorption signature against a bright 
background radio source, minihaloes will produce absorp- 
tion or emission against the CMB as well. The measured 
antenna temperature through a halo will be Tb{0) — 
T'cmb(O) exp(— Ti), where n is again given by Eq. (I13|l . where 
now Tcmb(2) = Tcmb{0){1 + z) is used for Tb{z) in Eq. Q, 
to give 



VCMB- (21) 



In principle, systems in absorption against the CMB 
would produce a 21cm forest. From Eq. (|A4|) for the post- 
shock temperature of a halo, however, it follows that re- 
quiring Tah < Tcmb{zc) restricts the range of halo masses 
in absorption to Af ^ 73OOAf0. This is comparable to the 
Jeans mass, so that the signal would be very weak. 

Similarly, the haloes could produce discrete emission 
lines against the CMB. In the limit Ts » Tcmb{zc), 
Eq. (pT]) becomes 

1/3 



15.07 



0.665 



M 



{nmh^y/s Vio6Af0 

/ M 



kHz 



1/3 



kHz, 



VTo^ ; ^^^^ 

independent of the redshift of the halo. The detection of 
an individual emission line in a narrow frequency channel 
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against a signal as weak as the CMB, however, is unlikely for 
the forseeable future. More viable is the average signal from 
an ensemble of haloes, integrated over several channels. The 
observed frame equivalent width weighted comoving cross- 
section of the haloes is given by 



~ 135.5^ 



M 



22.6 



M 



(lO^Mg 



© 

1/3 



kHz — kpc 



kHz — kpc 



(23) 



In the limit Ts >> Tcmb{zc) and using the 

Press-Schechter form for the halo mass function 

l|Pre88 fc Schechteij Il974l '). the optical depth Eq. (fTS)) 
may be cast in the suggestive form 



-—A ^tH(zc) T« 

327r H{Zc) rCMB(2:c) 



X [erf (in 



erf(tniin)] : 



(24) 



where T(Af) = (1 — Y)Qt/^m was assumed. The maxi- 
mum and minimum values oft, however, must be adjusted to 
match the halo mass function inferred from numerical simu- 
lations. The corresponding observed brightness temperature 
differential is 



5Tb{0) 



(4.6mK)n-'/2(l + Ze) 



sl/2 



(25) 



1 



-1/2 



where fm = [erf (tmax) ^erf (tmin)] is the mass fraction of the 
Universe in minihaloes. The signal from minihaloes is thus a 
small fraction of the available signal from the IGM. At red- 
shifts z < 20, when coUisional decoupling from the CMB is 
weak within the diffuse IGM, minihaloes may dominate the 
signal until Lyman resonance line radiation becomes avail- 
able to sufficiently decouple the spin state of the hydrogen 
atoms in the diffuse IGM from the CMB. The signal from 
the diffuse IGM, however, will dominate over that of the 
minihaloes once the mass-averaged 21cm efficiency of the 
diffuse component exceeds the mass fraction of the Universe 
in minihaloes. 

As shown by the expressions above for the number of 
absorption features along a line of sight and their effective 
optical depth, the ensemble statistics of the 21cm features 
depend on the number den sity of minihaloes. The fitting 
formula of lReed et al.l (12007') is adopted throughout the re- 
mainder of this paper. This choice is justified in Appendix iBl 



3 THE DYNAMICAL MINIHALO MODEL 
3.1 Model parameterization 

Spherical minihalo models for intergalactic gas clouds with 
the gas confin ed by dark matte r haloes were introduced by 
llkeuchil l|l986l ) and iReed (| 19861 ). with the lower mass end 
limited by that required to bind photoionised gas. Recogniz- 
ing that cosmological haloes are in gene ral dynamica l , with 
the gas either collapsing or expanding, iBond et al.l (|l988l ') 
extended the models by building the haloes from grow- 
ing spherical cosmological density perturbations. Since the 
haloes form from the merger of smaller haloes, there is some 
freedom in the choice of the initial linear gas density profile. 



Whilst a tophat perturbation is simplest, iBond et all (|l988l ) 
chose a Gaussian profile as it was found to best represent the 
final virialised halo density profile. In this paper, a tophat 
initial density profile is adopted for its simplicity, since the 
gas density will respond to the gravitational potential of the 
dark matter and gas combined, which is not very sensitive 
to the underlying density profile of the dark matter. 

The temperature and residual ionisation of the back- 
grou nd unper turbed IGM are solved for using RECFAST 
( Sea ger et al. I [2000). Fits to the temperature and residual 
ionization fraction are provided in Appendix [XJ along with 
estimates for the expected post-collapse minihalo tempera- 
tures. 

As galaxies or active nuclei form in the Universe, the 
radiation from these systems may heat the IGM tempera- 
ture to higher values. The hydrodynamical models are run 
without this heat input to provide a reference point for min- 
imal heating. The effect of a warm diffuse IGM on the 21cm 
signatures, however, is a lso ex plored. 

The approach of B ond et al.i f|1988 ) is used to describe 
the collapse of spherical mass shells, with the equations of 
hydrodynamics solv ed for numerica lly. In this paper, the nu- 
merical methods of lMeiks"inl l| 19941 ) are used. The gas tem- 
perature is solved for including atomic coUisional and radia- 
tive processes, as described in Mcijksin (1994), and cooling 
by molecular hydrogen. The ionization state of the gas is 
computed as well in order to account for hydrogen recombi- 
nations within the halo cores. The computation of molecular 
hydrogen formation and cooling is described in Appendix [Cl 

At the low mass end, the baryonic fluctuations are re- 
stricted by thermal pressure support, providing a lower limit 
to the minihaloes which contribute to the 21cm signal be- 
fore the fluctuations dissolve into sound waves. It is shown 
below that the high mass end is limited by the onset of star 
formation. 



3.2 Lower limiting mass 

The fluctuations in the baryonic component are filtered at 
short length scales by the homogenising effect of sound waves 
on scales comparable to the Jeans length. For linear per- 
turbations, the Jeans length is Xj = Cs{n /Cpm)^^^ , where 
Cs is the sound speed in the gas and pm is the total mass 
density (|Peebleslll980l ). Using the kinetic temperature from 
Eq. (|A1|) . this corresponds to a proper Jeans length and total 
halo mass Mj = (47r/3)pM(2:)(Aj/2)'^ for adiabatic fluctua- 
tions at redshift z of 



Aj ~ 2.2(H-2)-°-^2^ kpc; 



M, 



220(1-^2) 



Mq,(26) 



for 6 < 2 < 60. 

The Jeans mass is only an approximation to the lower 
limiting mass. Lower mass haloes will still contribute a non- 
negligible 21cm absorption signal against a bright back- 
ground radio source distinguishable from the diffuse IGM. 
As discussed in Sec. 12.21 above, the cumulative absorption 
signal is dominated by the lowest mass haloes still able to 
retain their gas. 

Haloes with masses near the Jeans mass will contribute 
an emission signal against the CMB provided Tsh > Tcmb, 
or M > 7300 M0, although with a reduced 21cm efficiency. 
An illustrative computation at 2 = 20 is shown in Fig- 
ure [1] for a 2.0 x IO'^A/q halo, comparable to the Jeans 
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Figure 1. Evolution of tlie fluid variables for a 2 X 10* Mq top 
hat spherical perturbation collapsing at Zc = 20. The Jeans mass 
at this epoch is Mj ~ 1.7 X 10* Mq. Shown are the hydrogen 
density (top left panel), gas temperature (top right panel), fluid 
velocity (bottom left panel) and 21cm efficiency ?JCMB (bottom 
right panel), all as functions of comoving radius. The curves corre- 
spond to 2 = 50 (dashed line), z = 30 (short-dashed long-dashed 
line), 2 = 20 (solid line) and z = 15 (dotted-dashed line). The 
peak temperature at z = 20 matches the predicted binding tem- 
perature Tbintj ~ 119 K (see text). At the epoch of collapse, 
VCMS < at r > 0.8 kpc, corresponding to absorption against 
the CMB rather than emission. 



mass Mj ~ 1.7 x 10* M0. The 21cm efficiency is appreciable 
within the core. In fact, the halo mass is below the critical 
mass required for a shock to form, as given by Eq. (|A7|) . The 
gas temperature is determined instead by the adiabatic flow 
of the gas as it falls into the halo and establishes hydrostatic 
equilibrium, settling at the binding temperature, Eq. (|A5|) . 
within the core. At z = 20, the temperature is sufficiently 
low for r > 0.8 kpc (comoving) that, from the 21cm effi- 
ciency factor 77CMB given by Eq. (|16p . the halo will absorb 
relative to the CMB rather than emit at these radii. The 
mass-weighted efficiency within the virial radius, however, 
still corresponds to net emission, although much reduced 
from the core value, with (jycMs) — 0.15. Because of the 
large number of such small haloes, their net contribution to 
the 21cm signatures is non-negligible, particularly at high 
redshifts. 



3.3 Upper limiting mass 

It has long been suggested that the first collapsing 
Jeans unstable gaseous structures in a Big Bang cosmol- 
ogy wou l d ha v e masses on the orde r of 10* — 10^ Mq 
l|Gamowl 1 19481 : [Peebles fc Dickd 1 19681 ). with star forma- 
tion triggered by molecular hydrogen cooling following 
the formation of H2 through gas phase processes with 
H~ and Hj" acting a s catalysts ( | Saslaw fc Zipovl 1 19671 : 
[Peebles fc Dick3 Il968l : iHirasawal Il969l ). This picture has 



survived remarkably well in contemporary cosmological 
models dominated by cold dark matter, with compara- 
bly small masses for the first collapsing objects inferred 
jBond fc Szalavll 19831: [Blumenthal et al.lll984l:|Peebleslll984 



Blurnenthal et al.ll 19851: lOstriker fc Gnedinlll996l:lAbel et all 



2OOOI: [Fuller fc CouchmanI I2OOOI : iMachacek et al.l l20o'ir 
Reed et al.ll2005l ). 



To account for star formation, the creation of molecu- 
lar hydrogen through gas phase reactions is computed dur- 
ing the collapse (see Appendix [C} . Only the formation via 
H~ is included, as t he added molecular hydrogen forme d 
via H+ is negligible (|Palla et al.lll983l : iLepp fc Shulllll984l ). 
When the cooling time is shorter than the c haracteristic in- 
flow time, the gas will be the rmally unstable l|Malagoli et al.l 
[l987; Balbus & So kei1ll989l ). Overdense pockets of gas are 
then assumed to cool rapidly a nd are removed from the 
flow isochorically l|Meiksinl 1 19881 ). When the cooling time 
is shorter than the inflow time, mass is removed from the 
flow and converted into stars at the rate 



(27) 



where tcooi is the net cooling rate within a gas shell, and 
is a dimensionless efficiency coefficient of order unity. 

Once an adequate mass of stars accumulates within a 
halo, it is presumed that a sufficient number of massive stars 
will have formed to photoionise the cloud or drive a wind 
through it via the mechanical energy input from wind losses 
and supernovae and completely disperse the gas in the cloud. 
For a Salpeter stellar initial mass function with a minimum 
stellar mass of IMq and maximum mass of lOOM©, an in- 
stantaneous burst with a total stellar mass of lO^M© formed 
will produce an ionizing luminosity of 10*^'* phs~^ for a du- 
rati on of 3.5 x 10^ yr, a ssuming a metallicity of 0.05 so- 
lar jLeitherer et al.lll999l ). Sufficient photons would be pro- 
duced to photoionize 6x10® Mq of hydrogen, and likely lead 
to the photo-evaporation of the gas in the halo. 

Even if radiative recombinations radiated away almost 
all the photoionization energy, typically at least one star as 
massive as 50Mq will also be produced, corresponding to a 
massive O-star. Such a star has a lifetime of under 10® yr, af- 
ter which it will explode as a Type II supernova with a char- 
acteristic mechanical energy input of 10*^ erg. Spreading 
the energy over a characteristic halo baryon mass of lO'' Mq 
corresponds to a gas temperature of Tsn — 7 x 10* K. Even 
allowing for 90% of the energy to radiate away, this is more 
than adequate to unbind the gas. A threshold of IG^Mq is 
therefore used as a criterion for sufficient star formation to 
disperse the gaseous content of the halo. It is found that 
this quantity of stars typically forms soon after the halo 
collapses. It is also found that the epoch and minimal halo 
mass for forming sufficient stars to disrupt the halo is not 
very sensitive to the assumed parameters, such as the star 
formation efficiency q, or the minimum required threshold 
of stars formed. 

More generally, an upper limit to the halo mass be- 
low which star formation will disrupt a halo may be es- 
timated as follows. After 10^ yrs, the formation of M« — 
10'' Mq of stars with masses between 1 — 100 Mq will 
produce approximately lO*'^ erg of mechanical energy in 
the form of winds and supernovae ejecta (ranging from 
]^q51.7 _ ]^q5i.9 g j ,^ £qj. q gg gQjg^j- gglar metal abundances) 

l|Leitherer et al.lll999l ). Allowing for a fraction /heat of this 
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energy to go into heating the gas, and for a fraction of 
baryons for a total halo mass M, setting (3/2)fcBT'hoat/7Ti = 
(/heat//6)(£^3/M)(M./lOOOM0), where S3 is the mechani- 
cal energy input from the formation of 10^ Mq of stars, gives 
a characteristic gas temperature to which the baryons will 
be heated of 

Theat ^ 2.95 X 10VhcatS3,52M.,3M6-' i^^^^ ^' ^^^^ 

where S3, 52 = Ss/lO^^ erg, M.,3 = M,/10^Mq and Mg = 
Mq. Comparison with the binding temperature from 
Eq. (jASP gives for the minimum halo mass required to retain 
the baryons (Theat < Tbind), 

M6,^i„ = {J^^l El%Ml';\l + z.)-''\{29) 

For /heat = 0.5, at Zc = 20 this requires a minimum halo 
mass of M > 1.5 x lO'^ Mq, corresponding to a post-shock 
temperature of T^h ^ 9200 K from Eq. HA4p . This mass is 
well above the halo masses found below required to form 
10^ Mq of stars. 

Previous estimates of the 21cm signature have not in- 
cluded the role played by molecular hydrogen cooling, pre- 
suming that molecular hydrogen is dissociated by a meta- 
galactic UV radiation field produced by early stars, so that 
star formation become s dominated only by more massive 
(10'^ - 10* Mq) haloes (|Haiman et al.lll997al ). Since it is the 
lower mass haloes that form first, and the gas within them 
does not survive long once the first stars form within them, 
it is unclear whether a sufficiently strong radition field will 
develop at early epochs. Nonetheless, both models with and 
without molecular hydrogen formation are considered. 

An illustrative computation including molecular hydro- 
gen formation is presented in Figure[2]for a 1.3 x 10® Mq halo 
collapsing at 2 = 20. The corresponding post-shock temper- 
ature is 1770 K from Eq. HA4|) . well below the temperature 
required for coUisional ionisation. The central temperature 
in the halo lies below this value because of efficient molec- 
ular hydrogen cooling. A comparison run with molecular 
hydrogen cooling turned off gives a central temperature of 
T ~ 2040 K. The evolution of the H2 fraction for the run 
is shown in Figure O along with the stars formed following 
H2 cooling. The column density is computed along lines of 
sight through a sphere defined by the turn-around radius at 
each epoch, except at 2 = 50, for which the comoving ra- 
dius of the initial perturbation is used. It is found that by 
2 = 20.05, shortly before the collapse epoch, W^Mq of stars 
has formed for g« = 1. 

The entropy generation by the shock just within the 
turnaround radius results in an inverted entropy profile, as 
shown in Figure [2] (The entropy per particle is computed as 
s = kBlogiT'-^^^/nn).) The inflowing gas may become con- 
vectively unstable and the gas turbulent within the core, 
although the timescale Iby = [{2/3){g/kB)\ds/dr\]~^^^ (the 
Brunt- Vaisala timescale), where g is the gravitational accel- 
eration, is long at these radii. At 2 = 20.05, it is shortest 
at r 10 kpc (comoving), where tav — 69 Myr. After this 
time, a spherically-symmetric computation will no longer 
be able to follow the evolution of the system in detail. Be- 
cause the mechanical energy input by the most massive stars 
formed will eject the gas on a shorter timescale, however, the 
halo will be evacuated before becoming convective. Haloes 
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Figure 2. Evolution of the fluid variables for a 1.3 X IO^Mq 
top hat spherical perturbation collapsing at Zc = 20. Shown are 
the hydrogen density (top left panel), gas temperature (top right 
panel) , fluid velocity (bottom left panel) and entropy per particle 
(in units of fe^ ) (bottom right panel) , all as functions of comoving 
radius. The curves correspond to 2 = 50 (dashed lino), z = 30 
(short-dashed long-dashed line), and z = 20.05 (solid line), by 
which 10^ Mq of stars has formed. The peak temperature at 2 = 
20.05 lies somewhat below the predicted post-shock temperature 
of ~ 1770 K because of efficient molecular hydrogen cooling. 



that reach rapid cooling prior to their collapse redshifts will 
contribute negligibly to the 21cm signature by the time the 
dark matter collapses. This forms a natural upper limit to 
the mass range of the haloes that do contribute. 

The central temperatures found for the collapsing 
haloes are shown in Figure |3] for a range of halo masses and 
redshifts. At low masses, the values with and without molec- 
ular hydrogen cooling both agree well with the expected 
post-shock temperature, as given by Eq. (|A4[) . Whilst the 
temperature without molecular hydrogen cooling continues 
to agree well with the expected post-shock temperature at 
higher masses, though tending towards the virial tempera- 
ture, the temperature reaches a ceiling of T ^ 1400 K when 
molecular hydrogen cooling is included. 

For central temperatures T ^ 8500 K, atomic cool- 
ing through the coUisional excitation of Lya becomes 
efficient, and stars are presumed to form according to 
Eq. (|27p . as for the case with molecular hydrogen cooling. 
The breaking of the scaling relation with mass predicted 
by Eq. (|A4[) results in a departure from cosmological self- 
similar accretion expected i n an Einstein-deSitt er universe 
(IFillmore fc GoldreichI [l98i : iBertschinged Il985h . The rea- 
sons for the departure are examined in § 15.61 

In the presence of molecular hydrogen cooling, the max- 
imum halo mass for which fewer than 10'^ AIq of stars form 
prior to the collapse of the halo is fit to 10% accuracy over 
8 < 2c < 30 by 
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Figure 3. Evolution of tiie H2 fraction (top left panel) and den- 
sity of stars formed (bottom left panel) for a 1.3 X 10^ Mq top 
hat spherical perturbation collapsing at Zc = 20, as a function 
of comoving radius. The curves correspond to 2 = 50 (dashed 
line), 2 = 30 (short-dashed long-dashed line), and z = 20.05 
(solid line), by wrhich IO'^A/q of stars has formed. The mass- 
weighted mean 21cm emission efficiency within the virial radius 
is ('?cmb) — 0.90 at the final time. During its growing phase, 
the perturbation would appear in absorption against the CMB 
(bottom right panel). 



Figure 4. Central temperature as a function of halo mass. Re- 
sults from numerical computation are shown at 2 = 10 (squares), 
20 (triangles) and 30 (circles). Results including H2 cooling are 
shown as solid symbols; open symbols show results without H2 
cooling. The curves show the expected post-shock temperature 
from Eq. | |A4| |. at z = 10 (short-dashed line), 20 (long-dashed 
line) and 30 (solid line). Molecular hydrogen cooling restricts the 
core temperature to T < 1400 K. With molecular hydrogen for- 
mation suppressed, the central temperature of increasingly more 
massive haloes approaches the virial temperature. 



Miooo ^ 10' f t4^) Mq. (30) 

This agrees well with fully three-dimensional hydrodynam- 
ical simulations of the minimum halo mass of typically 
10^ — 10^ Maj required for star f ormation at high redshifts 
dAbel et al. 2000; Fuller fc Couc hman 2000; Machacck ct al. 
I2OOII ). 

If molecular hydrogen formation is suppressed, higher 
mass haloes in which Lya cooling becomes efficient are re- 
quired for stars to form. In the absence of molecular hydro- 
gen formation, the maximum halo mass for which fewer than 
10^ Mq of stars form prior to the collapse of the halo is fit 
to 10% accuracy over 8 < < 50 by 

M^SJS^ ~ 9.1 X 10** exp[-(l + 2,)/51] Mq. (31) 

These masses are about a factor 5 greater than found when 
molecular hydrogen cooling is included. 

3.4 21cm signature against a bright radio source 

The absorption against a bright background source (BBS) is 
computed from Eq. (|4]). An illustration is given in Figure [5] 
for a halo with mass M = 0.9 x 10® Mq, including II2 cool- 
ing. (The results for the case with II2 cooling suppressed are 
nearly identical for this mass.) The absorption features are 
found not to evolve rapidly for a given halo mass, showing 
only a mild increase in the line-centre optical depth with 
time. The increase results from the continual inflow of ma- 



terial and the departure from self-similar accretion, as dis- 
cussed below. The trend of increasing optical depth with 
decreasing collapse epoch is contrary to the redshift scaling 
expectation for an idealised tophat halo, Eq. (|18|l , which pre- 
dicts a decreasing optical depth. The inflow also produces 
a substantially broader and shallower feature than occurs 
when the velocity broadening is suppressed. 

The profiles shown in Figure [5] were generated by inte- 
grating only along distances within the virial radius of the 
minihalo. An observed feature would include the absorption 
by the IGM in the surrounding gas. In Figure [6] profiles are 
shown at several impact parameters integrating out to the 
diffuse IGM. The intergalactic diffuse component is clearly 
visible in the wings of the features beyond a frequency offset 
of 4 kHz (in the observed frame). Since individual features 
would be measured relative to the diffuse absorption level, 
the offset at 4 kHz serves as a useful operational definition 
for the baseline relative to which the equivalent width of the 
feature would be determined. This definition for the absorp- 
tion equivalent width against a bright background source 
will be used throughout the remainder of this paper. 

According to Eq. (|lip . the contribution of a halo of mass 
M to absorption features with an observed equivalent width 
exceeding a given w'^'^ scales inversely with the mean free 
path Ainfp('ii'°o'') for intercepting a halo along a line of sight 
with observed equivalent width exceeding w"^" . The mean 
free paths are shown for haloes with H2 cooling collapsing 
at 2;c = 8 in Figure [7] for equivalent width limits wSq'^ > 0.1, 
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Figure 5. Evolution of the absorption feature against a BBS 
along a line of sight through the centre of a halo of mass M = 
0.9 X 10® Mq, for collapse epochs of = 15 (short-dashed line), 
10 (long-dashed line), and 8 (solid line). Profiles are shown both 
including the broadening from the line-of-sight peculiar veloci- 
ties (heavy lines), and without the peculiar velocity contribution 
(light lines). The frequency offset from line centre is in the ob- 
served frame. 



Figure 7. Comoving mean free path for haloes collapsing at Zc = 
8 (solid symbols) and Zc = 10 (open symbols) for cross sections 
giving rise to absorption features against a BBS with observed 
equivalent widths w°^^'^ > 0.1 kHz (triangles), to°Jj'= > 0.2 kHz 
(squares) and w'^'^ > 0.4 kHz (circles). 

0.2 and 0.4 kHz. Most of the systems with wZ]^" > 0.1 arise 
from haloes in the mass range 4 x 10'* — 3 x 10^ Mq. 
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Figure 6. Optical depth as a function of observed 21cm frequency 
for a halo of mass M = 0.9 X 10® Mq and collapse epoch of z^ = 8. 
Results shown for lines of sight at a) 6x = Oi b) bx = 1 kpc (co- 
moving), c) 6x = ^t.a. and d) bx = 2rt.a., where n.a. is the 
turnaround radius of the halo. The optical depth is computed 
through a sphere extending to the linear regime of the perturba- 
tion. Shown are the signal through the sphere (solid lines), and 
after subtracting off the values at 4 kHz (dashed lines). 



3.5 21cm signature against the CMB 

The collapsing haloes will produce emission features rela- 
tive to the CMB. Whilst the signal from an individual halo 
will not be detectable for the forseeable future, their profiles 
reveal the degree to which the signal of a single minihalo 
may be distinguished from that of the surrounding gas. The 
temperature differentials in the observed frame are shown in 
Figure |8] for a halo with mass M = 0.9 X 10® Mq, including 
H2 cooling. (The case with H2 cooling suppressed differs lit- 
tle for this mass.) The increase in the brightness temperature 
for later collapse epochs arises from the continual inflow of 
material. As in the case for absorption against a bright back- 
ground source, the infall broadens the feature and reduces 
its amplitude compared with the resuff when the velocity 
contribution is suppressed. 

Emission against the CMB is not very sensitive to the 
signature from the infalling gas surrounding the minihalo. 
This is because collisions are able to decouple the hydro- 
gen spin structure from the CMB only at sufficiently large 
overdensities near a collapsed halo. The equivalent width 
weighted cross section E^*"^ is only moderately affected by 
extending the integration volume beyond the core radius, 
as shown in Figure [O] Since the mean cosmic optical depth 
scales in proportion to this factor, as in Eq. (|13p . the sur- 
rounding gas is not expected to contribute much more to the 
mean cosmic optical depth. It is the truncation of the sig- 
nal with impact parameter that provides the justification for 
treating the minihalo contribution to the total IGM 21cm 
signature as distinct from the diffuse component. 
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Figure 8. Evolution of the observed brightness temperature rela- 
tive to the CMB along the line of sight through the centre of a halo 
of mass M = 0.9 x 10*' Mq\, for collapse epochs of Zc = 15 (short- 
dashed line), 10 (long-dashed line), and 8 (solid line). Profiles are 
shown both including the broadening from the line-of-sight pe- 
culiar velocities (heavy lines), and without the peculiar velocity 
contribution (light lines). The frequency offset from line centre is 
in the observed frame. 



Figure 10. Comoving equivalent width weighted cross section (in 

kHz-kpc^) for emission against the CMB, as a function of halo 
mass for haloes collapsing at Zc = 10 (squares) and z^ = 20 (cir- 
cles). The cross section is shown in the observed frame. Shown for 
models including H2 cooling (solid symbols) and without (open 
symbols). Also shown are fits for models without H2 cooling at 
Zc = 10 (solid line) and Zc = 20 (dashed line). 
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Figure 9. Observed integrated equivalent width weighted co- 
moving cross section for emission against the CMB as a function 
of comoving maximum impact parameter 6™^", normalized by 
the core radius. Shown for haloes of mass M = 0.9 X 10*' Mq 
(solid hues), M = 0.3 X 10*' Mq (long-dashed lines), and M = 
0.8 X 10^ Mq (short-dashed lines), for collapse epochs of Zc = 8 
(heavy lines) and 15 (light lines). The turnaround radius corre- 
sponds to b^'^^'/rc ^ 3.89. 



Figure 11. Ratio of the equivalent width weighted cross section 
for emission against the CMB for the hydro dynamical models to 

the cross sections for the tophat model, as a function of halo 
mass for haloes collapsing at Zc = 10 (squares) and Zc = 20 
(circles). Shown for models including H2 cooling (solid symbols) 
and without (open symbols). 
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The geometrically-averaged observed emission equiva- 
lent width is shown for a range of halo masses for haloes 
collapsing at Zc = 10 and 20 in Figure [TO] The steep rise in 
the mean equivalent width with halo mass is partly a con- 
sequence of the mass scaling given by Eq. (|22p . but is also 
due to the mass dependence of the 21cm radiative efficiency 
?7CMB of the haloes. The more massive a halo, the greater 
its 21cm radiation efficiency for emission against the CMB. 

Polynomial fits to the observer-frame equivalent- 
width weighted halo cross sections of the form EjJ'^ = 
Ef^o'JiClogio Af)*"' are shown as well. The fits are per- 
formed for halo masses logjQ M > 5, below which E^''^ is 
vanishingly small. Table [T] provides the coefficients for a 
range of collapse epochs. The fits may be used to interpolate 
the results of the hydrodynamical models on both the halo 
mass and the collapse epoch. 

The ratios of the equivalent width from the hydrody- 
namical model to the tophat model for the haloes are shown 
in Figure 1111 The scaling of the equivalent width with halo 
mass approaches the tophat prediction for the larger mass 
haloes. The average values from the hydrodynamical models 
are much smaller than the tophat prediction from Eq. (|22p 
at the low mass end, where the haloes approach the Jeans 
mass and the gas becomes less overdense, as shown in Fig- 
ure[T] The hydrodynamical values agree well with the tophat 
prediction at the high mass end despite the very difi'erent 
internal structures of the haloes. The smallest mass haloes 
actually produce a net absorption signal (cf. Figure [1} . 



4 COSMOLOGICAL 21CM SIGNATURE 
STATISTICS 

4.1 Detection of minihaloes 

The models neglect 21cm absorption or emission from the 
large-scale surroundings of the minihaloes. As primordial 
perturbations are intrinsically ellipsoidal (BBKS), and flat- 
ten as they grow (^Peebles,, 1980i' ) , modelling the region around 
a minihalo is not straightforward. Only in the late stage of 
collapse will a spheroidal system form. Nonetheless, it may 
be expected the infall region surrounding the m inihalo may 
still be reasonably spherical (|Zhang et al.lll99^ ). 

The degree of spectral isolation of a minihalo may be 
estimated from the mean free path for intercepting a halo 
out to the turnaround radius. Approximating the space den- 
sity of the haloes as dn/d\og M gives a mean free path 
for interception of a typical halo of mass 0.9 X IO'^Mq, 
with a comoving turn- around radius of 12.5 kpc, collaps- 
ing at = 8 of 3.4 Mpc (comoving). The mean (observed) 
frequency separation of these systems is about 200 kHz. 
The detection of an individual absorption feature against 
a bright background source would require measuring the 
brightness temperature differential within frequency chan- 
nels not much greater than 8 kHz, and preferably as narrow 
as 4 kHz or smaller to resolve the feature. For a spectral sig- 
nal smoothed over broader frequency intervals comparable 
to the frequency separation of the features, the systems will 
produce a small modulation of the intensity step from the 
diffuse gas. For smoothing over much wider frequency in- 
tervals, it follows from Eq. p3|l that the contribution from 
all systems will simply add a small constant amount to the 
overall intensity step from the diffuse gas. 
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Figure 12. Top panel: Minimal Lya photon collision rate, com- 
pared with the thermalization rate Pthi required for absorption 
against the CMB of the diffuse component of the IGM to exceed 
that from minihaloes, allowing for H2 formation in the haloes 
(solid line) or not (dashed line). Bottom panel: The 21cm effi- 
ciency rjcMB for absorption or emission against the CMB from 
the diffuse component of the IGM due to H-H atomic collisions 
alone. For no heating (solid line), the IGM is in absorption, and 
dominates the minihalo contribution for z > 19. For an IGM 
temperature = 2Tcmb (2) (short-dashed line) and IOTqmb (^) 
(long-dashed line), the diffuse IGM is in emission, and dominates 
the emission from minihaloes for 2 > 15. Also shown is the mini- 
halo mass fraction Jm of the IGM, allowing for H2 formation in 
the haloes (heavy dotted line), or not (light dotted line). 



It is noteworthy that at large redshifts absorption by gas 
beyond the turnaround radius diminishes T:"^^ somewhat for 
emission against the CMB, especially for small mass haloes. 
As shown in Figure [9] the asymptotic value of E°'''^ has not 
been reached for a M = 0.8 x 10^ Mq halo collapsing at Zc — 
15, so that the temperature differential may be somewhat 
over-estimated at these redshifts. 

When the 21cm efficiency of the diffuse component of 
the IGM exceeds the mass fraction of the IGM in minihaloes, 
the diffuse IGM signal against the CMB will overwhelm the 
minihalo signal. Given the mass range of minihaloes that 
contributes to the minihalo signal, it is possible to estimate 
the critical Lya scattering rate required for the diffuse IGM 
signal to dominate. Using Eq. (|16|l . and assuming a minimal 
halo mass given by M^h and a maximum mass of Miooo, or 
Af "q(^^ if H2 formation is suppressed, produces the critical 
scattering rates shown in Figure [12] (upper panel) . Only a 
small fraction of the thermalization rate is required. In fact, 
for 2 > 19 collisions alone are sufficient to ensure absorption 
from the diffuse IGM dominates the signal from minihaloes 
(lower panel) . In the presence of moderate amounts of heat- 
ing, the diffuse IGM will dominate over minihaloes at even 
smaller redshifts. Shocks in the diff use IGM may produce 
further emission against the CMB (iFurlanetto fc Loebll200j : 
iKuhlen et al.ll2006l : IShapiro et al.ll20oi r 
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Table 1. Coefficients of S°''=(M) 
kHz-kpc'^ (comoving). 



S*_gai(logjg A/)* ' for emission against the CMB, for models witliout H2 cooling. Units are 
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Figure 13. Observed equivalent width distributions for absorp- 
tion against a BBS. Top panel: Shown are the differential distri- 
butions at 2 = 8 (solid lines) and z = 10 (dashed lines), with 
H2 cooling (heavy lines) and without (light lines). The dotted 
lines show fits through the distributions with H2 cooling. Bottom 
panel: The corresponding cumulative distributions. 



4.2 21cm statistics for absorption against a bright 
radio source 

The most readily measured statistic of the minihalo 21cm 
signature against a bright background source is the equiv- 
alent width distribution. The distributions at z = 8 and 
z = 10 from the hydrodynamical computations are shown 
in Figure 1131 The cases with and without H2 cooling have 
very similar distributions, except that for z — 8, the higher 
mass haloes that survive with H2 formation give rise to a 
slightly higher frequency of high equivalent width systems, 
with w°^^" > 0.32 kHz, although these systems are rare. 

The distributions are well-fit in the range 0.1 < 



A^'o exp(- 



< 0.3 kHz by the exponential curve N / dvo'^^dz 



7w.). At 



8, A^o 



1400 and 



0.069 kHz, while at z = 10 iVo = 3100 and w. = 0.062 kHz, 
both for the case with H2 cooling. 

The distribution is dominated by very low mass haloes 
only for observed equivalent widths below w"^^ ^ 0.02 kHz. 



Figure 14. Observed brightness temperature differential relative 
to the CMB as a function of redshift, including H2 cooling (solid 
line) and without (dashed line). 



Excluding haloes with masses no greater than an order of 
magnitude above the Jeans mass severely suppresses the 
number of absorption features for such small equivalent 
width values. 



4.3 21cm statistics for emission against the CMB 

The evolution of the observed brightness temperature is 
shown in Figure [T?] Since the emission equivalent width in- 
creases with the mass of the haloes, as shown in Figure [TOl 
the extension to higher halo masses when H2 formation is 
suppressed produces a doubling of the differential tempera- 
ture compared with the case including H2 cooling. 

At high redshifts, much of the temperature differential 
arises from low mass haloes. Haloes with masses less than an 
order of magnitude above the Jeans mass contribute about 
one third of the signal at z = 20, both for the cases with 
and without H2 cooling. The contribution of the low mass 
haloes, however, is much reduced by z = 8. 

The cumulative observed emission equivalent width dis- 
tributions for haloes collapsing at Zc = 8 and 10 are shown 
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Figure 15. Observed equivalent widtli cumulative distribution 
for emission against the CMB at z = 8 (solid lines) and z = 15 
(dashed lines), with H2 cooling (heavy lines) and without (light 
lines). 

in Figure [15] The heavy curves show the distribution includ- 
ing H2 coohng and the light curves without. The extension 
to higher mass haloes in the latter case produces a higher 
frequency of larger equivalent width systems, although the 
distributions are found to extend to w°^^ > 5 kHz for both 
cases. Numerous absorption systems with —0.01 < w°^^ < 
kHz arise as well for Zc = 10 resulting from absorption in 
the outer regions of the accreting haloes. 



5 DISCUSSION 

5.1 21cm signature statistics 

The 21cm signatures of minihaloes considered here, absorp- 
tion against a bright background radio source and emis- 
sion against the CMB, are dominated by haloes in differ- 
ent mass ranges. The absorption signal against a bright 
source is dominated by haloes in the approximate mass 
range 4.5 < logjg(M/Af0) < 6. The less massive haloes 
in this range lie near the Jeans mass, so that the post-infall 
gas pressure prevents as large an overdensity within the halo 
core from developing as in the more massive haloes. Haloes 
with masses above this range have a reduced absorption ef- 
ficiency as a result of their higher post-shock temperatures. 
By contrast, the emission signature against the CMB in- 
creases with increasing halo mass, limited only by the max- 
imum halo mass before star formation leads to the ejection 
of the halo gas. 

As a consequence of their differing halo mass depen- 
dences, the signals will be affected differently by the various 
parameters defining the models. In this section, the role of 
some of these parameters in determining the strength of the 
signals is examined. 

An issue only partially addressed in the minihalo ap- 



proximation is the contribution of the surrounding moder- 
ately overdense gas on the overall signals. Minihaloes are 
not generally isolated structures, but are located within 
larger scale density inhomogeneities which themselves are 
non-linear. In the case of absorption against bright back- 
ground sources, the overdense gas surrounding the collapsed 
minihalo was shown to contribute a small amount to the sig- 
nal. Because of its near uniformity in frequency, it could be 
readily subtracted off, leaving the equivalent width largely 
unaffected. 

In the case of emission against the CMB, the role played 
by the surrounding gas is less clear. Beyond the central emit- 
ting minihalo, the gas is sufBciently dense and cool to act 
as an effective absorber against the CMB, cancelling in part 
the emission in an unresolved observation. If the absorbing 
region is sufficiently extensive, it could substantially reduce 
the net emiss i on fr om minihaloes, as has been argued by 
iKuhlen et all l|2006l) on the basis of numerical simulations. 
It is found here that extending the line-of-sight integration 
around the minihalo, although producing some reduction in 
the overall equivalent width, does not much reduce it. It is 
important to ensure the simulations adequately resolve the 
minihaloes in extended structures that would appear to be 
absorbing if under-resolved numerically. Achieving the re- 
quired resolution in a simulation volume sufficiently large 
to represent a fair sample of the universe, without itself be- 
ing overdense, is a formidable computational challenge. A 
possible solution is to model the statistical fluctuations sur- 
rounding the minihaloes semi-analytically. 

5.2 Sensitivity to cosmological parameters 

A prominent factor in establishing the strength of the sig- 
nals is the primordial power spectrum. Both 21cm signa- 
tures are sensitive to the amount of power on small scales 
through their dependence on the halo mass distribution. In- 
deed, the signatures probe scales smaller than any current 
measurement of the primordial power spectrum, providing 
a potentially powerful means of constraining the shape of 
the primordial power spectrum and even the nature of dark 
matter. 

In Figure [TUl the predicted cumulative observed equiv- 
alent width distributions for absorption against a bright 
background source are shown for the best-fitting cosmo- 
logical models including a running spectral index found by 
combining WMAP Year 7 CMB data with the data from 
the A tacama Cosmology Telescope (ACT) jPunklev et al.l 
I2OIOI ). Two models are considered. The first is defined by 
the cosmological parameters Qm ~ 0.330, fit, = 0.670, 
Qbh^ = 0.02167, h = 0.661, ag^-i = 0.841, spectral in- 
dex at fco = 0.002 Mpc~^ of n{ko) = 1.032, and slope 
dn/d\ogk = —0.034. The second incorporates additional 
statistical priors based on measurements of the Hubble con- 
stant and Baryonic Acoustic Oscillations (BAOs), and is 
given by the cosmological parameters f2,„ = 0.287, ^Iv = 
0.713, Qth'^ = 0.02206, h = 0.691, erg,,-! = 0.820, spec- 
tral index at fco = 0.002 Mpc~^ of n(fco) ~ 1.017, and slope 
dn/d log fe = -0.024. 

The suppression of power on small scales severely re- 
duces the expected number of systems, by as much as two 
orders of magnitude for the first model. Including the con- 
straints from Ho and BAO measurements substantially in- 
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Figure 16. Cumulative observed equivalent width distributions 
for absorption against a BBS, at 2; = 8 (solid lines) and z = 10 
(dashed lines), with H2 cooling (heavy lines) and without (light 
lines). Top panel: Cumulative distribution predicted for the best- 
fitting running spectral index model constrained by ACT and 
WMAP7 data. The dotted lines show the fit distributions with H2 
cooling for the fiducial cosmological model. Bottom panel: Cumu- 
lative distribution predicted for the best-fitting running spectral 
index model constrained by ACT, WMAP7, Ho and BAG data. 
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Figure 17. Observed brightness temperature differential relative 
to the CMB as a function of redshift, for cosmological models 
including a running spectral index. Shown are curves using models 
constrained by ACT and WMAP data with additional constraints 
from measurements of Hg and BAOs (solid lines), and without 
(dashed lines), both for minihalo models including H2 cooling 
(heavy lines) and without (light lines). 



Figure 18. Observed central equivalent width for emission 
against the CMB as a function of halo mass, at 2 = 10 (squares), 
15 (triangles) and 20 (circles). Results including H2 cooling are 
shown as solid symbols; open symbols show results without H2 



cooling. The dashed lines show power laws w° 



creases the expected number of systems, but the numbers 
still lie well below the predictions for the fiducial cosmo- 
logical model used in this paper. The counts of absorption 
systems would readily distinguish between the two running 
spectral index models, illustrating the power the 21cm de- 
tection of minihaloes has for constraining the small-scale 
primordial density power spectrum in the absence of other 
suppression mechanisms of the absorbers. 

The reduced small scale power also reduces the pre- 
dicted temperature differentials from the CMB by factors of 
a few to several, as shown in Figure 1171 Comparison with 
Figure [14] shows the reduction is particularly pronounced at 
high redshifts. Measurements of the temperature differen- 
tial especially at high redshift could thus readily constrain 
the running spectral index in the absence of other sources 
of suppression of the signal. It is noteworthy that a very 
weak signal at high redshifts compared with that expected 
in the fiducial model could be mistaken for a detection of cos- 
mic reionisation in the absence of other observational con- 
straints. 

The halo mass dependence of the equivalent width for 
emission against the CMB given by Eq. (|22|l suggests the 
possibility of measuring the halo mass function directly 
from 21cm measurements. Whilst the hydrodynamical mod- 
els produce a diminishing equivalent width with impact pa- 
rameter, an observation with sufficient angular resolution to 
resolve the core of the structure could establish the central 
equivalent width value. The central value is shown as a func- 
tion of halo mass in Figure 1181 The equivalent width varies 



with the halo mass approximately as Wi, 



This is 



steeper than predicted for the tophat model. Nonetheless, 
the trend is sufficiently tight to yield a useful measurement 
of the halo mass. 
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The opposing mass dependences of the equivalent width 
as measured in absorption against a bright background 
source (BBS), Eq. p9p . and the emission equivalent width 
in the tophat model suggests combining them may eliminate 
the dependence on halo mass, and yield a direct measure- 
ment of the combination of cosmological parameters given 
by 



(obs \ / obs \ 



4.42(gi^(kHz)^ (32) 



The ratio of equivalent widths eliminates the dependence on 
the baryon fraction, and would yield a measurement of the 
halo mass through 



(■ obs\ 

v"^o J gym 

y^vn /BBS 



M 



\ 106 Mg 



2/3 



(33) 



In fact, dissipation breaks these relations. The actual re- 
lations both are mass dependent, and have no real advan- 
tage over inferring the mass from a measurement of the core 
equivalent width using the scaling relation shown in Fig- 
ure [181 Probing the halo cores in emission against the CMB 
would, of course, be a formidable observational challenge, re- 
quiring a space-based facility or possibly one placed on the 
far side of the Moon to achieve both the required sensitivity 
free of radio Earthshine (radio-frequency interference) and 
the long baseline. The larger minihalo cores would just be 
resolvable within the limitations imposed by image blurring 
by interplanetary and interstellar plasma, with a cha racter- 
istic rms spread of {{SOsfY''^ ~ l(i//lMHz)~2 deg (|Carillil 
I2OO8I ). 



5.3 Effects of first light 

The far ultra-violet light from the first radiation sources 
in the Universe will redshift into local Lya resonance line 
photons sufficiently distant from the sources. The scattering 
of the photons off the neutral hydrogen serves as a mecha- 
nism to decouple the spin temperature of the hydrogen from 
the CMB temperature through the Wouthuysen-Field effect 
(|Wouthuvsenlll952l : lFieldlll958h . Prior to reionisation, the 
intensity will grow sufficiently strong that the Lya colli- 
sion rate Pa will approach the thermalization rate Pth- For 
a cold IGM, with Tk < Tcmb, the 21cm optical depth will 
increase, increasing the absorption signature of minihaloes 
against bright background radio sources. 

Well before the epoch of reionization, a small scatter- 
ing rate Pa « Pth is expected. Whilst the photons will 
only slightly decouple the spin temperature from the CMB 
temperature, the vast reservoir of neutral hydrogen outside 
collapsed haloes will produce a sizeable absorption signature 
for an IGM colder than the CMB (Figure [T2)) . 

Another possible effect of radiation from the first 
sources is to heat the IGM. Prior to reionisation, the 
primary heating mechanisms are through the photoelec- 
tric absorption of x-rays generated in shock-heated col- 
lapsed structures, within supernovae remnants from early 
evolved massive stars, or possibly from any x-ray stellar 
sources or Active G a lactic Nucle i (AGN) that have formed 
(Madau et al.' 1997; Tozzi ct al. 2000; C ariUi et all |2002| : 
[Purlanetto fc Loeb 2004; K uhlcn c t al. 200^, or by the scat- 
tering of higher order Lyman resonance line photons in the 



vicinity of their sources (*Mei ksinll2010l ). The amount of x-ray 
heating is unknown. Because both 21cm absorption against 
a bright radio source and the signal against the CMB are 
sensitive to the temperature of the IGM, they provide a 
means of probing the early heating of the IGM by the first 
radiation sources. Several estimates have been made in the 
literatur e for the evolutio n of the IGM temperature. The 
model of lFurlanett3 (|2006l ) is representative, with a heating 
rate based on the expected amount of star-formation at high 
redshifts, scaling by the low redshift correlation between x- 
ray luminosity and star-formation in starburst galaxies. The 
diffuse IGM temperature is found to cross the CMB temper- 
ature in the redshift range 12 < z < 15 for plausible model 
parameters. By z = 10, the IGM temperature is expected to 
exceed 100 K, and may exceed 1000 K. 

An increase in the IGM temperature will have two ef- 
fects on the 21cm signatures: an increase in the minimum 
halo mass giving rise to a signal and the strength of the sig- 
nal from individual haloes itself. The degree of the hydro- 
dynamical effect on the haloes will depend on the specifics 
of the heating rate and history. In particular, as discussed 
in Appendix the minimum halo mass required for the 
gas falling into the halo to shock and virialise increases. In 
lower mass haloes, the gas is heated by adiabatic compres- 
sion, and the halo gas density profiles will adjust accordingly. 
The fraction of haloes affected will increase gradually with 
the increase in the IGM temperature. 

Heating the IGM to a temperature above that of the 
CMB will produce a qualitative difference in both the ab- 
sorption signature against a bright background radio source 
and emission against the CMB since for both cases, dis- 
lodging the spin temperature from the CMB temperature 
through either collisions or Lya photon scattering will in- 
crease the spin temperature. To isolate this latter effect, the 
spin temperature is re-computed assuming an instantaneous 
boost of the IGM temperature by ATk ~ IOTcmb , and the 
impact on the 21cm signatures assessed. A set of hydrody- 
namical models with gradual heating is also computed to 
assess the impact of the hydrodynamical response of the gas 
on the signal. 

5.3.1 21cm signature against a bright radio source 

The spin temperature for the hydrodynamical models is re- 
computed for five cases, having Pa/Pth ~ 0.001, 0.01, 0.1, 
1 and 10. Figure [19] shows typical profiles through a halo of 
mass M = 0.9 x 10^ Mq collapsing at Zc = 10, including H2 
cooling, at a (comoving) projected separation from the halo 
centre of b± = 0.9 kpc, corresponding to an observed equiv- 
alent width = 0.20 kHz for Pa = 0. For Pa = O.lPth, 
the equivalent width is only slightly reduced. For Pa — Pth, 
the equivalent width is nearly halved, to 0.11 kHz. The pro- 
file is distorted into one having two minima, suggestive of 
two overlapping features. This arises from an enhanced con- 
tribution from the outflowing gas beyond the turn-around 
radius, as shown in Figure 1201 Most of the absorption at 
V — = 1 kHz arises from infalling gas near the centre of 
the halo {xios < 12.5 kpc), however the Doppler boosted 
absorption by outflowing gas also contributes. Whilst the 
contribution from the outflowing gas is small for Pa = 0, for 
which the spin temperature is coupled to the CMB temper- 
ature, for Pa = Pth the spin temperature is coupled to the 
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Figure 19. Absorption line profiles through a 0.9 X 10 Mq halo 
collapsing at Zc = 10 at a projected separation from the halo 
centre of bj^ =0.9 kpc (comoving), corresponding to an observed 
equivalent width ui°^^ = 0.20 kHz for Pa = 0. The profiles are 
for Pa/ P%ii = (solid lines), 0.1 (long-dashed lines) and 1 (short 
dashed lines). The profiles relative to the background level at 
V — VQ = A kHz are shown as light curves. The frequency offset is 
in the observed frame. 



Figure 21. Absorption line profiles through a 0.9 X 10 Mq halo 
collapsing at Zc = 10 at a projected separation from the halo 
centre of = 5.9 kpc (comoving), corresponding to an observed 
equivalent width w°^'' = 0.030 kHz for Pa = 0. The profiles 
are for Pa/Pn^ = (solid lines), 1 (long-dashed lines) and 10 
(short dashed lines). The profiles relative to the background level 
a,t V — Vf) = i kHz are shown as light curves. For Pa > fthi the 
profiles would appear as emission features against the background 
absorption level. The frequency offset is in the observed frame. 
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Figure 20. The differential optical depth dr/dl per unit comov- 
ing length along the line of sight (in comoving kpc) through a 
0.9 X 10® Mq halo collapsing at Zc = 10 at a projected separation 
from the halo centre of fej^ = 0.9 kpc (comoving), shown for fre- 
quency offsets V — UQ = (solid lines) and 1 kHz (dashed lines) in 
the observed frame. Top panel: Pa = 0. Bottom panel: Pa = Ptl\- 



much colder IGM temperature, enhancing the net absorp- 
tion to beyond that at line centre. 

At still larger projected separations, referencing the 
background continuum to the value at 4 kHz results in posi- 
tive apparent "emission" features for sufficiently large Pa , as 
shown in Figure [2T] A line of sight through the halo above at 



5.9 kpc corresponds to w° 



0.030 kHz for Pa 



0. 



For Pa = Pth, the absorption feature would appear as an 



emission line with w?, 



0.050 kHz relative to the back- 



ground absorption level at u - 
the feature would have w?}'^ : 



uo = 4: kHz. For Pa = lOPth, 
-1 0.12 kHz. Such strong mock 
emission lines against a bright background source would be 
a tell-tale signature of an intense Lya background radiation 
field and a cold IGM. 

For a Lya scattering rate as large as the thermalization 
rate, the optical depth increases substantially in the absence 
of heating, since the IGM temperature is much lower than 
that of the CMB. As a consequence, weak absorption fea- 
tures along the lines of sight adjacent to low mass haloes 
become substantially enhanced, resulting in a large num- 
ber of absorbers. The observed equivalent width cumulative 
distributions for haloes collapsing at Zc = 10 are shown in 
Figure [22] for Pa = Pth and Pa = lOPth (top panel). (The 
equivalent widths are computed using the intensity at a fre- 
quency offset of 2 kHz instead of 4 kHz as the continuum 
value due to limitations imposed by the volume of the hy- 
drodynamical computations for such high absorption optical 
depths when Pa > Pth-) The distributions are only represen- 
tative, as features will be produced by much more moderate 
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Figure 22. Cumulative observed equivalent width distributions 
for absorption against a BBS from haloes collapsing at Zc = 10, 
including H2 cooling. Top panel: The curves correspond to = 
Pth (solid line) and Pa = lOPth (dashed line). Bottom panel: The 
spin temperatures are computed including a temperature boost 
by 10rcMB(2c)- The curves correspond to Pa/Pth = (solid 
line), 1 (long-dashed line) and 10 (short-dashed line). Also shown 
(dot-dashed line) is the cumulative distribution allowing for the 
hydrodynamical response of the halo gas to a gradual heating of 
the IGM, without H2 cooling and for Pa = (see text). In both 
panels, the dotted line shows the fit distribution with H2 cooling 
for Pq = and no boost in the IGM temperature. 

density fluctuations throughout the IGM, not only in the 
vicinity of haloes. They indicate, however, the large num- 
bers of features that would arise in a cold IGM. 

The mean free path between haloes corresponding to 
absorption features with observed equivalent widths w'^'^ < 
-0.05 kHz and wZl" > 0.2 kHz and 0.4 kHz are shown in 
Figure [23] (top panel). Systems with w"^" > 0.2 kHz arise 
predominantly from low mass haloes, with a flat contribu- 
tion over the range 10* — 10^ h~^MQ. Mock emission sys- 
tems, with uijjjj'^ < —0.05, arise from the outer regions of 
more massive haloes with masses 6 - 9 X 10^ h-^Mco- 

In an IGM in which Pa was fluctuating, patches with Pa 
exceeding Pth would show up as regions with a highly dense 
21cm forest. Similarly, if Pa exceeded Pth throughout the 
IGM while reionization and pre-reionization heating were 
still incomplete, patches with low Tk would produce a highly 
dense 21cm forest. 

The effect of allowing for an instantaneous tempera- 
ture boost ATk ~ lOTcMB in the IGM on the absorption 
equivalent width distribution is shown in Figure [22l (bottom 
panel) for Pa — 0, Pth and lOPth- The effect is to severely re- 
duce the expected number of systems with w"]^" > 0.1 kHz. 
Weaker features arise predominantly from the peripheries of 
haloes in the mass range 2 — 5 x 10^ H^^AIq, as shown in 
Figure [221 (bottom panel). Colhsional coupling is weak for 
this gas, so that the spin temperature is strongly coupled to 
the CMB temperature and is little affected by a boost in the 



Figure 23. Comoving mean free path for haloes collapsing at 
Zc = 10 giving rise to absorption features against a BBS. Top 
panel: The mean free path for Pa = Pth (solid symbols) and 
Pa = lOPth (open symbols), corresponding to absorption sys- 
tems with observed equivalent widths w"^^ < —0.05 kHz (circles), 
«J°J^= > 0.2 kHz (triangles; upper limits shown as inverted trian- 
gles) and w°^^ > 0.4 kHz (squares). Bottom panel: The mean 
free path for haloes, adjusting the spin temperature by adding 
lOTcMB instantaneously to the kinetic temperature of the gas. 
Shown for Pa = (solid symbols) and Pa = Pth (open sym- 
bols), corresponding to absorption systems with observed equiv- 
alent widths w"^^ > 0.05 kHz (triangles; upper limits shown as 
inverted triangles), w°^^ > 0.1 kHz (squares) and w°^^ > 0.2 kHz 
(circles). 

IGM temperature. The dramatic steepening in the equiva- 
lent width distribution would be a means of identifying the 
presence of a temperature boost. 

The mean free path for a feature with w°^^ > 0.2 kHz 
increases from Amfp — GOh~^ Mpc (comoving) for ATk ~ 
and Ptj = 0, as shown in Figure [T] to Amfp — 200 Mpc 
(comoving) or larger for ATk ~ IOTcmb- The signal is dom- 
inated by haloes with a mass M ~ 3 X lO'^ /i"^Mq. Weaker 
systems, with w°^^ > 0.1 kHz, arise primarily from lower 
mass haloes, with 1 - 5 x IO'^/i'^Mq, for Pa = 0. For 
Pa — Pth, the signal is dominated by a nearly flat con- 
tribution from haloes with masses 2 x 10^-10^ h~^MQ. 

To assess the hydrodynamical impact of a gradual tem- 
perature boost, a set of models collapsing at Zc — 10, with- 
out H2 cooling, is computed including a constant photoelec- 
tric heating rate per volume by hard photons. The hard radi- 
ation field is turned on at 2 = 20 with an amplitude adjusted 
to increase the temperature of the diffuse IGM by IOTcmb 
at z — 10, corresponding to Tk — 302 K for the diffuse com- 
ponent. Adiabatic compression in a halo with a post-collapse 
density increase by a factor will heat the collapsed gas 
to /c^Tk — 9500 K. This greatly exceeds the post-shock 
temperature of T^h — 2100 K, from Eq. (|A4[l . for a halo as 
massive as Mh — 4 x 10® Mq. As a result, the gas temper- 
ature increases gradually within the core, and the central 
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Figure 24. Fluid variables for a 4 X 10® Mq tophat spherical 
perturbation collapsing at Zc = 10, without H2 cooling. Grad- 
ual photoelectric heating is added by a hard source turned on at 
2 = 20, with intensity adjusted to boost the diffuse IGM tem- 
perature by 300 K at z = 10. Shown at z = 50 (dashed line), 
15 (short-dashed long-dashed line), 10 (solid line) and 8 (dotted- 
dashed line). Shown are the normalised hydrogen density (top 
left panel), gas temperature (top right panel), fluid velocity ex- 
pressed as a Mach number (bottom left panel), and entropy per 
particle (in units of kg) (bottom right panel). An outflow from 
the core develops for z < 10. Also shown in all panels are the 
corresponding results for the same halo model at 2 = 10 without 
any photoelectric heating (dotted lines). 



density reaches a peak value of a factor of 50 smaller than 
for the case without external heating, as shown in Figure [24l 
The accretion does not establish a steady state, with an out- 
flow from the core developing by z = 8. The gas density is 
comparable to the dark matter density within the core at 
z = 10, so that a more realistic 3D computation may result 
in additional compression of the dark matter. The adiabatic 
heating term, however, would still play an important role in 
determining the hydrodynamical structure of the halo. The 
equivalent width along a line of sight passing through the 
centre of the halo is reduced by a factor of 2.6 compared with 
the model without external heating, and corresponds to an 
observed equivalent width of w"^^ ~ 0.13 kHz. Higher mass 
haloes computed are found to produce no larger equivalent 
widths. 



Figure 25. The observed brightness temperature relative to the 
CMB along the line of sight at a projected radius 6^ = Tc for 
a halo of mass M = 0.9 X 10® Mq collapsing at z^ = 15, for 
Pa/Pth = (solid line), 0.01 (short-dashed line), 0.1 (long dashed 
line), and 1 (dot-dashed line). The emission signature converts 
into a net absorption signature by Pa/Pth = 0.1. The frequency 
offset from line centre is in the observed frame. 
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5.3.2 21cm signature against the CMB 

The scattering of Lya photons weakens the 21cm emission 
signature from minihaloes against the CMB as the spin tem- 
perature couples increasingly strongly to the lower kinetic 
temperature of the difluse IGM. In Figure [25l the brightness 
temperature signature is shown for Pa/Pth ~ 0, 0.01, 0.1 and 
1. The corresponding observed equivalent widths are 9.8, 
4.1, -42 and -270 Hz, integrating between f — i/o = ±4 kHz. 



Figure 26. Observed integrated equivalent width weighted co- 
moving cross section for emission against the CMB as a function 
of comoving maximum impact parameter b^^^ , normalized by the 
core radius. Shown for a halo of mass M = 0.9 X 10® Mq collapsing 
at Zc = 15, for Pa/Pth = 0-001 (dotted line), 0.01 (short-dashed 
line), 0.1 (long-dashed line) and 1 (solid line). 
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The signal is halved for Pa/Pth ~ 0.01, and produces a net 
absorption signature for Pa/Pth ~ 0.1 and larger. 

In fact the equivalent width of an individual halo is no 
longer well-defined, since the entire IGM becomes absorbing. 
From Eq. (|13p , the strength of the overall signature depends 
on the equivalent width weighted cross-section Tt"^^. Except 
for very small values Pa/Pth << 1, even if the minihalo con- 
tribution were defined by restricting the frequency range for 
evaluating the equivalent width to ±4 kHz of the line cen- 
tre, T,'^^ is again not well-defined for an individual halo. As 
shown in Figure [26l averaging over increasingly large maxi- 
mum impact parameters produces an increasingly negative 
value even for Pa/Pth as small as 0.01. The emission signal 
from the minihalo is swamped by the absorption from the 
surrounding cold IGM. Indeed, on larger scales absorption 
from the diffuse IGM is expected to dominate, as shown in 
Figure 1121 For even small values of Pa/Pth, the minihalo 
model is no longer an effective approximation for estimating 
the signature of the IGM against the CMB. If a direct signal 
from the minihaloes may be detected at all, it would only be 
from high angular resolution measurements that were able to 
resolve individual minihaloes. This would require subarcsec- 
ond beam angles, well beyond the specifications of currently 
existing or planned radio facilities for the forseeable future. 

It is possible that the Lya scattering rate would be 
suppressed within dusty regions produced by winds from 
galaxies, as the Lya photons would be absorbed along their 
long path lengths as they re-scatter. In the vicinity of a 
bright galaxy, however, the scattering of higher order Lyman 
resonance line photons alone would be adequate to produce 
a net absorption signature from nearby minihaloes. Typical 
scattering rat es of P„/Pth — 0.001 — 0.01 for Ly-n photons 
are expected (|Meiksinll2O10l ) . 

Boosting the IGM temperature by ATk- = IOTcmb 
gives rise to an emission signature from the diffuse IGM, 
which produces a wing on the minihalo emission line, as 
shown in Figure 1271 A substantial emission wing is found 
even for Pa = Q. The coUisional coupling to the gas ki- 
netic temperature alone is sufficient to produce the wing 
(cf . Figure I12p . Even when the frequency range for evalu- 
ating the equivalent width is restricted to within 4 kHz of 
the line centre, the equivalent width weighted cross-section 
SSJ^" is found not to converge with increasing impact pa- 
rameter out to at least a few times the turnaround radius, 
as shown in Figure [28l Whilst there is a quasi-convergence 
for b± ~ rc, averaging over larger values results in a di- 
verging quantity. The emission signal from the surroundings 
of the halo swamps that of the minihalo contribution from 
fax < rc- The minihalo model is again no longer a useful 
approximation for quantifying the contribution to the total 
temperature differential Tb — Tomb , now dominated by the 
contribution from the diffuse IGM beyond the core radius 
of the minihalo. 
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Figure 27. The observed brightness temperature relative to the 
CMB along the line of sight at a projected radius 6^ = ''c for 
a halo of mass M = 0.9 X 10^ Mq collapsing at Zc = 15 and 
including a boost in the IGM temperature by ATff = IOTcmb i 
for Pa/Pth = (solid line), 0.01 (short-dashed line), 0.1 (long 
dashed line), and 1 (dot-dashed line). The growing wings arise 
from the increase in the diffuse IGM emission as Pa increases. 
The frequency offset from line centre is in the observed frame. 
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5.3.3 Combined 21cm signatures 

The impact of the first radiation sources on the joint 21cm 
signals from the IGM prior to reionisation is divided into 
six broad regimes: (a) ATa- << Tk, Pa <<< Pth; (b) 
ATk < TcMB - Tk, Pa «< Pth-, (c) ATk < Tcmb - Tk, 
O.OlPth < Pa « Pth; (d) ATk > Tcmb - Tk, Pa « Pth; 
(e) ATk « Tk, Pa > Pth; and (f) ATk > Tcmb - Tk, 



Figure 28. Observed integrated equivalent width weighted co- 
moving cross section for emission against the CMB as a function 
of comoving maximum impact parameter b™'"', normalized by 
the core radius. Shown for a halo of mass M = 0.9 X 10® Mq 
collapsing at Zc = 15 with the IGM temperature instantaneously 
boosted by ATk = IOTcmb, for Pa/Pth = (soUd line), 0.01 
(dotted line), 0.1 (short-dashed line) and 1 (long-dashed line). 
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Pa > Pth- The diflerences in behaviour between the absorp- 
tion against bright background radio sources and the signal 
against the CMB may be exploited to distinguish the fac- 
tors of heating, Lya scattering and the uncertain amount 
of power on small scales. 

In regime (a), absorption by minihaloes against a bright 
background radio source and emission against the CMB is 
unaffected by galactic feedback. (By Pa «< Pth is meant a 
Lya scattering rate well below that required to increase the 
21cm efficiency of the diffuse component of the IGM to above 
the mass fraction of the IGM in minihaloes. See Figure [121 ') 
The redshift evolution of the signals traces the growth of 
structures on comoving scales down to a few kiloparsecs. 
In regime (b), the absorption signal against a bright back- 
ground source weakens whilst the minihalo signal against the 
CMB is largely unaffected. Although the signal from shocked 
gas in the diffuse IGM may evolve somewhat, it may be pos- 
sible to isolate the development of large-scale structure using 
the signal against the CMB from the effects of heating on 
the absorption systems against bright sources. In regime (c), 
the signal from the absorption systems weakens as in regime 
(b), but the signal against the CMB now goes into deep ab- 
sorption. This would indicate the presence of a moderate 
metagalactic UV radiation field. Disentangling the effects 
of heating, Lya scattering and small-scale power would be 
difficult. In regime (d), the absorption against a bright back- 
ground source continues to diminish as the IGM is warmed, 
but the signal against the CMB goes into strong emission, 
which continues to strengthen as either ATk or Pa increases 
until the signal saturates with full emission from the IGM 
against the CMB. Before saturation, joint modelling of the 
signal against bright sources and against the CMB may par- 
tially disentangle the effects of heating, Lya scattering and 
small-scale structure growth as an increasing Pa has little ef- 
fect on the absorption sytems but a strong effect on the emis- 
sion signal against the CMB once Pa > O.OlPth- In regime 
(e), the absorption signal against background sources be- 
comes very strong, producing a large number of absorption 
systems, whilst the signal from the minihalos and their envi- 
rons against the CMB will go strongly into absorption. The 
large number of absorbers indicates little heating so that 
its strength is determined primarily by the Lya scattering 
rate and the amount of small-scale power. In regime (f), the 
absorption signal from the minihaloes against a bright back- 
ground source will be diminished, with a large contribution 
arising from the diffuse IGM. The diffuse IGM signal against 
the CMB will be in emission and saturate for a sufficiently 
large Lya scattering rate. The minihalo signal against the 
CMB will be overwhelmed by the emission from the diffuse 
component. 

5.4 Effects of reionisation 

The collapsed haloes achieve large central hydrogen den- 
sities and large central H i column densities, with values 
A'h > 10^° cm"'^ typical. The haloes will be detectable into 
the Epoch of Reionisation, although with reduced cross- 
sections as their outer regions become heated and pho- 
toionised by the ambient UV radiation field. The degree of 
reduction depends on the intensity of the radiation field. 

Two simple limiting cases are considered for estimating 
the incident fiux of ionising photons. The emissivity neces- 



sary for reionising the Universe may be minimally estimated 
by requiring one photon per hydrogen atom over a Hubble 
time. If clumping drives the mean hydrogen recombination 
time to under a Hubble time, then the emissivity may be 
estimated instead by requiring one photon per mean hydro- 
gen recombination time. The column density required for 
shielding the interior of the halo from the incident fiux may 
then be estimated by modelling the ionisation zone within 
the halo as an inverted Stromgren sphere. 

The minimal fiux based on one photon per baryon is 
{nu{z))c, where the angle brackets indicate spatial averag- 
ing over the IGM. If the radiation field penetrates to a depth 
/ within the halo, balancing recombinations with ionisations 
within the ionised surface layer gives / ~ {nu{z))c/n^^^aB, 
where nn.c is the internal hydrogen density of the halo and 
a_B the radiative recombination rate (Case B) within the 
ionised zone. A density of nu.c ~ fc^inu) and a temper- 
ature of 10'* K in the ionised layer gives A'h = nn,cl ~ 



X 10^" cm"^ for = l/(187r 



2x1/3 



Since the 

density rapidly rises within the collapsed halo to higher val- 
ues, the required column density will be somewhat smaller, 
depending on the mass of the halo and collapse epoch. 

If the dumpiness of the IGM is sufficiently large, the 
reionisation will be recombination limited. In this case the 
incident flux of ionisation radiation is (n^ ^asc/ H{z). If 
the clumping factor of the IGM is dominated by the haloes, 
then balancing recombinations within the outer ionisation 
layer of a halo with the ionisation rate gives a depth for the 
ionisation layer of I = fvc/H{z), where fv is the volume 
filling factor of the haloes. For collapsed tophat haloes, this 
is related to the mass fraction /m in collapsed haloes by 
fv = fcfM- The column density of hydrogen through the 
ionised layer is then Ah ~ fMnB_{z)c/H{z) ~ lO'^^ cm~'^ at 
z = 10, adopting /a/ from Figure fT^ This means the ionizing 
radiation would be able to penetrate deeply into the core, 
much reducing the cross section, but still not eliminating all 
the haloes from detection. Many reionization scenari os have 
been considered in the literature IC ncdin 2000; Ciardi et al-l 
I2OO3I : iMellema et al]|2006l : IZahn et al.ll2007l ). 



5.5 Illustrative detections by radio facilities 

The most straightforward means of detecting the 21cm sig- 
nature of minihaloes is through absorption against a bright 
background radio source. Two detection strategies are pos- 
sible, the direct detection of individual absorption lines and 
the detection of the ensemble fiuctuations produced by the 
systems. For an absorption feature extending over Ach fre- 
quency channels of width Ach, the rms error on the equiva- 
lent width is 



1/2 



A, 



1/2, 



(34) 



where cta^i, is the rms instrumental noise level per channel 
and Jc is the continuum level of the background radio source. 

The expected variance in the measured transmissivity 
Q(A) — J dX' R{X — A') exp(— r>j/), where R describes the 
response function of the telescope, is given by 



{{SQf) = {{Q~{Q)f) 



(Qf 
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in the hmit that the response function is broad compared 
with the absorption features (jPress et aLiil993i '). Here, H{X) 
is given by 



r oo 

H{\) = / dro< 
X [ dX' [l 



dwx dro dX 



where d^JV/dwxdrodX is the number density of absorp- 
tion systems per line-centre optical depth tq per equiv- 
alent width (in wavelength) wx per rest wavelength A, 
X = {X — Ao)/AAd, where AXo — (6/c)Aio, and (j>xix) — 
■K~^^^ exp(— a;'^) is the dimensionless line profile. For ro << 
1, this simplifies to H{X) ~ / dTodwx{d'^ M / dwxdTodX)w\. 
Denoting the number of absorption systems per equiva- 
lent width per unit rest wavelength by d^M /dwxdX = 
[{l + z)/XQ]d^M/dwxdz, Eq. JaS]) becomes (still in the limit 
ro << 1), 



. 1 



1^0 



A 



(37) 



where 

poo , , 

Jo d."^' 



Wx 



XoWiyg/i'o was used and 



ur. ^ „h., iivV,"" ) was defined. The instrument re- 

spouse has been approximated here as a square-well of 
width A: _R(A) = 1 for \X\ < A/2 and otherwise. The 
rms fluctuation in Q due to instrument noise is ag = 

(Ach/A)l/VA,J/c = rTA/Ic 

The rms fluctuation {{SQYY^'^ over a band of width A 
arising from the minihaloes may be distinguished from the 
rms noise fluctuations ctq by comparing [{((5Q)^) +o"q1^''^ — 
aq with (tq/(2A^s)^^^, where gq is estimated by sampling 
the variance over A'^s bands of width A, corresponding to 
an rms on the rms of (jq / {2NaY^'^ . In this way, even when 
{{5QY) is small compared with cjq, the fluctuations due to 
the minihaloes may be detected at a significant level. For a 
significance level ua, the measured fluctuations are related 
to the rms fluctuation in the equivalent width measured with 
channels of width Ach for the same integration time by 



Its- obs-i2\l/2 
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where Am = A'^chAch is the width of the absorption feature 
and At = A^sA is the width of the band over which the noise 
level is sampled. 

The number of bright sources at z > 6 is unknown. 
Crudely extrapolating the comoving number density of the 
highest redshift known radio sources predicts ~ 1 — 2 x 10* 
sources per unit redshift with /c > 6 mJy over 8 < z < 15, or 
at least a few dozen out to 2; = 10 assuming a steep decline 
with redshift similar to that of bright quasars (|Carilli et al.l 
12002! ) . The SKA will have sufficient spectral resolution and 
sensitivity to detect individual absorption lines in such 
sources. A SKA pathfinder like LOFAR could do so as well 
for plausible source brightnesses at high redshifts, but the 
required integration times would be substantially longer. 

For a 24 hr integration, the SKA is expected to achieve 
an rms noise level of ai ~ 34 /ijy in a 10 kHz channel 
for frequencies in the range 100 — 200 MHz (|Carilli et al.l 
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Figure 29. The rms fluctuations arising from absorption by mini- 
haloes against a 6 mJy background radio source at 2 = 10, includ- 
ing the instrumental noise contribution, for SKA assuming 10 kHz 
channels and a 24-hr integration time. Shown is the case without 
feedback, allowing for molecular hydrogen formation (heavy line; 
solid points) and without (light line; open points). Also shown 
are several cases with feedback for the model with molecular hy- 
drogen formation at 2 = 10: an instantaneous boost in the IGM 
temperature by AT = IOTcmb with Pa = (solid square) and 
Pa = Pth (solid triangle). The star shows half the value for the 
case with Pa = Pth a^nd AT = 0. Also shown is the case for a 
gradual increase in the IGM temperature of AT = IOTcmb by 
2 = 10, including the hydrodynamical response of the gas for 
a model without H2 formation, with Pa = (open square) and 
Pa = Pth (open triangle). The straight solid line and dashed lines 
correspond, respectively, to the rms noise level of the telescope in 
10 kHz width channels and the rms fluctuations about the noise 
level averaged over 10"^ samples in a lOMHz wide band. 



|2002| ). The noise rms in the measured transmissivity will be 
1/"* (JQ = 0.0057 for a 6 mJy source. The fluctuations in the 
(38)transmissivity due to the minihaloes in a 10 kHz band is 
found to be ((^Q)^)^/^ ~ 0.0023 from Eq. §7} and using 
the equivalent width distribution of Figure [13] ai z — 8 al- 
lowing for H2 cooling, and {{SQf)'^^^ ~ 0.0028 without. A 
region As — 10 MHz long contains Ns — 10'^ samples of 
10 kHz width channels, so that ag may be measured to an 
accuracy of aQ/{2N^y^^ ~ 0.00013. The onset of the fluc- 
tuations may then be detected at a signiflcance level of 3.5a 
for the case with H2 cooling, and 5.2a without. 

For a feature extending over A'ch = 4 channels, Eq. ((34} 
(or Eq. (EHl)), gives {{Sw°^''fy^^ ~ 0.036 kHz, so that 5a de- 
tections may be made of absorption systems with observed 
equivalent widths w^^^ > 0.18 kHz. Thus the integration 
time required to detect the rms fluctuations in the trans- 
missivity arising from the minihaloes is sufficient to detect 
the stronger absorption features individually. For a longer 
integration time of 300 hrs, it would be possible to detect 
systems down to w°^^ > 0.05 kHz, enabling essentially the 
full equivalent width distribution to be measured. 
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In fact, even for the shorter integration time of 24 hrs 
the effect of the weaker absorption systems will have been 
detected as well. Over 90% of the contribution to {{SQ)'^)^^^ 
arises from systems with w°^'' < 0.1 kHz. A comparison of 
{{dQ)'^)y^^ with the number count of absorption systems 
with w°^'^ > 0.18 kHz would provide an estimate of the 
steepness of the equivalent width distribution, and so con- 
strain the amount of feedback to the IGM. As shown in 
Figure 1291 the fluctuations are not expected to evolve very 
rapidly out to z = 10 in the absence of feedback. Feedback, 
however, may substantially alter the evolution. Allowing for 
an ambient Lya radiation field producing a scattering rate 
Pa = Pth, with a negligible heat input, will dramatically 
boost {{5Q)^)^^^ , using the equivalent distribution of Fig- 
ure [H 

By contrast, in the absence of a Lya radiation field, 
an instantaneous boost in the IGM temperature by AT = 
lOTcMB will only marginally reduce {{5Q)^)^^^, as shown 
in Figure 1291 despite the large decrease in the number of 
systems with w"^^^ > 0.1 kHz. This is because the systems 
with < 0.1 kHz are little affected by the temperature 
boost since they arise from gas strongly coupled to the CMB 
temperature. Allowing for the hydrodynamical response of 
the gas for a gradual increase in the IGM temperature of 
AT = lOTcMB by 2 = 10 (see Sec. [SXTjl . shown for the 
case with no H2 cooling, similarly results in only a moder- 
ate suppression of the signal. Increasing Pa to Pth weakens 
both signals as it brings the spin temperature of the gas 
everywhere into closer equilibrium with the gas kinetic tem- 
perature. 

The much higher instrumental noise level of LOFAR 
will require a much longer integration time to detect the fluc- 
tuations, however only a fraction of the full 48 MHz band- 
width available would be required, allowing the remainder 
to be used for alternative projects, such as a deep sky sur- 
vey. As an illustration, for a noise level of a a = 6.0 mjy at 
150 MHz in a 12 kHz band for a 1-hr integration timeQ a 
5a detection of fluctuations of the magnitude {{5QYY^^ — 
0.0028 against a 6 mJy source at z = 8 would require an 
integration time of 1300 days. Should a source as bright as 
30 mJy be found, the integration time shortens to a more 
manageable 1200 hrs. For this integration time, it would be 
possible to detect individual absorption systems four 1 kHz 
channels wide with w?}^'^ > 0.18 kHz at the 5cr level. 



5.6 Departure from self-similar accretion 

The characteristic infall velocity and temperature of the 
post-shock halo gas predicted from Eqs. (IA3P and (|A4|I (Ap- 
pendix|3, scale as power laws of the halo mass. Self-similar 
accretion may then be expected, as is the case for cosmologi- 
cal i nfall of a self-gravita ting initially pressureless coUisional 
gas (|Bertschingerl [l985h . A departure from self-similarity, 
however, is produced by the outer boundary condition of a 
finite IGM temperature, and so a non-vanishing pressure. 
Whilst for more massive haloes the IGM pressure may be 
neglected, the pressure results in a non-negligible contribu- 
tion to the post-shock temperature in a minihalo. 



Based on data available from lofar.strw.lcidcnuniv.nl 
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Figure 30. Fluid variables for haloes collapsing at Zc = 20 for 
halo masses logjQ M = 6.44 (dotted line), 6.54 (dashed line), 6.63 
(dot-dashed line), 6.76 (short-dashed long-dashed line), and 6.83 
(solid line), for an initial tophat spherical perturbation, with H2 
cooling suppressed. Shown are the normalised hydrogen density 
(top left panel), gas temperature (top right panel), fluid velocity 
expressed as a Mach number (bottom left panel), and entropy per 
particle (in units of kg) (bottom right panel). The peak value of 
T(r)/Tsh and the entropy increase with halo mass. 



The jump condition for the temperature for an adi- 
abatic shock is T2/T1 = [2-yMl - (7 - 1)][(7 - 1)MI + 
2]/[(7-|- l)'^A^i] for a gas with ratio of specific heats 7, where 
Ml = |'!'i|/ci is the Mach number for the infiowing gas with 
velocity vi (relative to the shock front) and adiabatic sound 
speed ci. For an initially pressureless gas, the shock is al- 
ways very strong and the post-shock temperature depends 
only on the gas velocity, as in Eq. (|A4|I . For gas starting 
at the intergalactic temperature, adiabatic compression will 
pre- warm the gas by a factor of up to /^^, as given by 
Eq. (|A6|I . This results in a mild shock rather than a strong 
shock, with the strength of the shock increasing with the 
halo mass, as shown in Figure [5(71 For the larger masses, the 
post-shock temperature approaches the strong shock limit 
Tsh of Eq. HA4|I at about half the (proper) collapse radius 
rc = /c7-o/(l + Zc). 

As shown in Figure [30l the gas becomes self-gravitating 
as the gas continues to fiow inward. For the case of a halo 
with mass 6 x 10** Mq (solid curve), the baryonic mass 
reaches 7 x 10^ Mq in the inner 0.8 kpc (comoving) by 
z = 20.1, about 50 times the dark matter mass. As the 
gas comes to rest at the centre of the halo, the gas devel- 
ops a steep central density profile. The further compression 
in the time-varying potential results in a further boost in 
the gas temperature, until values somewhat above Tsh are 
reached at the centre. In more massive haloes, temperatures 
approaching, and in some cases even somewhat exceeding, 
the virial temperature are reached, driving the gas into a 
rapid cooling phase through Lya excitation, which will ul- 
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5.7 Comparison with previous results 
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Figure 31. Brunt- Vaisiila growth time for convectivc instability 
for haloes collapsing at Zc = 20 at 2 = 20. Shown for haloes 
of mass logjQ M = 4.73 (short-dashed line), 5.20 (dottcd-dashed 
line), 5.64 (long-dashed line) and 6.00 (solid line), as a function 
of radius, normalised by the core radius rc(M). 



timately result in star formation. More than 10'' Mq of stars 
form within the core by z = 20. The resulting supernovae 
will likely expel the gas from the minihalo. A cosmological 
simulation code unable to resolve the inner few hundred par- 
sees (comoving) would fail to detect the cooling instability 
that develops in the centre of the halo, and treat the halo 
as stable. 

One caveat in the solution is that the inverted entropy 
profile may result in a convective instability outside the core. 
As discussed in Sec. 13.31 above, the growth time is typi- 
cally longer than the lifetime of massive stars, so that the 
halo will be disrupted before the instability sets in. This 
will not be the case, however, for haloes of too small mass 
to form stars. As shown in Figure 1311 the Brunt- Vaisala 
timescale is shorter than 100 Myr just outside the core ra- 
dius for the more massive haloes less massive than the criti- 
cal mass logj^Q M = 6.1 for which more than 10^ Mq of stars 
form from H2 cooling. For the case log^ M = 4.73, close to 
the Jeans mass, the gas becomes convectively unstable just 
within the core. If the timescale is shorter than the typical 
time between mergers, the halo would become convectively 
unstable. In this case, the halo gas would fragment, possibly 
bringing on a contraction of the halo, diminishing the num- 
ber of systems detected per unit redshift. A high resolution 
3D simulation would be required to assess the effect of the 
instability. 

A second caveat is that the baryons will draw in the 
dark matter into the central region, an effect not included 
here. As a consequence, the hydrodynamical flow of the gas 
must be included to describe the dark matter halo profile 
on sub-kiloparsec (comoving) scales. The increase in the 
amount of dark matter will strengthen the gravitational field 
in the centre of the halo, and enhance the cooling instability. 



Fluctuations in the 21cm signature are an expected conse- 
quence of the large-scale structure of the dark matter and 
baryo ns in the Universe (|Hogan &: Reed Il979l : iTozzi et al.l 
|2000D . The 21cm forest from small-scale cosmological struc- 
tur es was extrac t ed fro m numerical reionization simulations 
by ICarilli et all (|2002h . The IGM was heated by the re- 
ionizing radiation to a volume-averaged temperature well 
above that of the CMB, and had a baryonic mass per par- 
ticle of at best logjQ M = 5.7. The masses of the resolvable 
dark matter haloes would be over two orders of magnitude 
larger, exceeding the upper mass limit before stars would 
form. Most of the features in the simulated spectra in fact 
arose from moderately overdense filaments, the gas within 
which will not have been well-resolved into the minihaloes 
that give rise to the 21cm features. The simulations were 
also carried out in boxes too small to provide fair samples of 
the universe. The simulations nonetheless clearly established 
the viability of detecting the 21cm forest in the spectra of 
bright background radio sources. 

To overcome so me o f the limitations of simulations, 
iFurlanetto fc Loebl (|2002l ) modelled the minihaloes as 
spheres in hydrostatic equilibrium within a dark matter 
halo with a static isothermal core at the halo virial tem- 
perature and an infalling outer region at the IGM tempera- 
ture extending to the turn-around radius, beyond which the 
sphere parameters were adjusted t o match onto the d iffuse 
IGM. The IGM temperature from ICarilli et all l|2002l ) was 
adopted; at a = 10 the temperature is 1000 K. A model 
with an IGM temperature a factor ten lower was also con- 
sidered. In both cases, the IGM temperature is above the 
CMB temperature. No models with an IGM temperature 
lower than that of the CMB were considered. They allow for 
minihaloes up to a mass corresponding to a virial tempera- 
ture of 10^ K. Lya photon scattering rates of Pa/Pth — 1 
and 10 at 2 = 10 were adopted, as well as Pa = 0. They 
find dN/dz ~ 10 for w"^" > 0.1 kHz at 2 = 10 for then 
colder model with Pa — Pth- The model parameters cor- 
respond closely to those considered here with an instanta- 
neous IGM temperature boost and Pa ~ Pth- For the model 
here, dN / dz{w'^'' > 0.1 kHz) ~ 1, about an order of mag- 
nitude smaller. This may in part be due to the somewhat 
higher IGM temperature of 300 K for the boosted temper- 
ature model here, compa red with 100 K for the model of 
IFurlanetto fc Loebl (|2002l ). Their hotter model, with an IGM 
temperature of 1000 K, has dN/dz{w°]^'' > 0.1 kHz) ~ 1. In 
the model here, however, the line density declines faster for 
systems with > 0.3 kHz, for which dN/dz ~ 0.002; 

IFurlanetto &: Loebl l|2002t ) obtain 0.02, an order of magni- 
tude larger. The difference may arise from the lower trun- 
cation in the upper minihalo mass imposed here, allowing 
for star formation following H2 formation. Given the differ- 
ences between the models, the overall agreement is reason- 
ably good. Of course, as discussed in Section r5.3.1l allowing 
for the hydrodynamical response of the gas to gradual pho- 
toelectric heating corresponding to a temperature boost of 
AT — lOTcMB by 2 = 10 eliminates all absorption features 
with w°;;'= > 0.13 kHz. 

Using a static, non-singular, truncated isothermal 
sphere model with a temperature given by the virial tem- 
perature, and immerical simulations to estimate the halo 
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number density. Ishapiro et al. IHooel) compute the expected 
21cm emission from minihaloes assuming radiative feedback 
effects from the first light sources are negligible. They al- 
low for halo masses having virial temperatures up to 10'' K. 
At z = 8, they obtain a brightness temperature differential 
5Tb ~ Tb — TcMB — 3.5 mK at 2: = 8. Adjusting to their 
cosmological model parameters (similar to those assumed 
here, except the primordial power spectrum is untilted and 
normalized to (Tg,j-i = 0.9), the tophat model for the same 
mass range predicts 5Tb — 3.0 mK, very close to their value. 
As shown in Figure [TT] however, the tophat model generally 
exceeds the prediction from the dynamical minihaloes com- 
puted here. The dynamical halo model, with H2 formation 
suppressed, predicts the smaller value of 5Tb — 1.4 mK, 
allowing for masses up until star formation occurs as given 
by Eq. pip . Using the to phat model to extend to the up- 
per mass limit adopted bv lShapiro et al.] (|2006l ) boosts this 
value to 1.9 mK. Given the differences between the models, 
this is reasonably good agreement. The weaker value found 
here, however , only slightly exceeds the value of 1.5 mK 
IShapiro et aD (|2006h find for the emission from the diffuse 
IGM alone. Allowing for the lower upper mass limits ob- 
tained here suggests the emission from minihaloes is at most 
comparable to that of the diffuse IGM. 



6 SUMMARY AND CONCLUSIONS 

A spherical collapse model is used to characterize the 21cm 
absorption signatures of minihaloes against bright back- 
ground radio sources and the Cosmic Microwave Back- 
ground prior to the Epoch of Reionization. The model 
evolves an initially linear spherical perturbation by simulta- 
neously solving for the evolution of the dark matter using a 
shell code and the gas using a hydrodynamics code. Two sets 
of models, with and without molecular hydrogen formation, 
are computed. The resulting models self-consistently include 
inflow and an accretion shock, and match onto cosmological 
boundary conditions on large scales. Atomic and molecular 
radiative processes are included to compute the temperature 
structure of the haloes. A mass sink is added scaling like the 
local gas density and inversely with the cooling rate to avoid 
a thermal cooling catastrophe from developing in the halo 
cores in the presence of strong cooling. The mass removed 
is presumed to form stars. 

A maximum minihalo mass for giving rise to a 21cm 
signal is set based on the formation of an adequate mass in 
stars to expel the gas from a minihalo either through photo- 
evaporation or a supernova-driven wind. When molecular 
hydrogen is allowed to form, the maximum mass at Zc = 10 
is about 1.5 x IQp Mq, declining gently with the collapse 
epoch as (1 + Zc)~^^'^. When molecular hydrogen formation 
is suppressed, so that cooling is due to atomic processes 
alone, the ma:ximum mass rises to 7 x lO'^ A/q at z = 10, 
declining exponentially with redshift as exp[— (1 -I- 2c)/51]. 
The central temperature of the maximum mass halo ai z — 
10 is 6500 K, well below the virial temperature of lO** K 
often assumed. The temperature is adequate for initiating 
the cooling required to form sufficient stars to result in the 
expulsion of the gas from the halo. 

An inverted entropy profile near the formation of an 
accretion shock may produce a bouyancy instability in the 



infalling gas. The timescales are generally long, but could be 
as short as 50-100 Myr. If the gas becomes unstable to frag- 
mentation, the details of the gas temperature and molecular 
hydrogen formation would be modified. 

Estimates in the literature for the minimum UV radia- 
tion field required to suppress molecular hydrogen formation 
were based on the formation of molecular hydrogen within 
collapsed haloes that initially had none. It is found here 
that haloes collapsing as early as Zc = 30 would form molec- 
ular hydrogen cores that were optically thick to dissociat- 
ing radiation. As these systems merge into larger haloes, it 
is possible a substantial reservoir of self-shielded molecular 
hydrogen could collect in the cores of subsequent genera- 
tions of haloes. In this case, the amount of star formation 
suppression could be less than previous estimates. 

The 21cm statistical signatures are computed using a 
halo mass function normalized to high resolution numeri- 
cal simulations. An estimate for the mass function of haloes 
below the numerical resolution is made by extrapolating a 
model based on peak statistics. Matching the two estimates 
for resolved haloes, the extrapolations to lower mass haloes 
are found to agree to better than 15%. The statistics are 
computed using the cosmological model constraints from 
WMAP. Using cosmological constraints that take into ac- 
count the ACT CMB data as well results in a substantial 
decrease in the strength of the signals, suggesting the 21cm 
signature from minihaloes is a sensitive probe of the primor- 
dial density power spectrum on small scales. 

Two scenarios are explored, one with no large-scale 
galactic feedback on the IGM and one with. In the case 
without feedback, the IGM temperature and residual ion- 
ization following the recombination epoch are adopted for 
the initial linear perturbations. No further energy or exter- 
nal radiation field are added. 

A very large number of absorption systems is produced 
against a bright background radio source. At 2 = 8 — 10, 
about 10 systems per unit redshift with an observed equiva- 
lent width w°^^ exceeding 0.1 kHz are predicted for the fidu- 
cial cosmological model considered. The numbers decline to 
a few per unit redshift for observed equivalent widths ex- 
ceeding 0.25 kHz. Renormalizing to the running spectral in- 
dex models constrained by both WMAP and ACT data de- 
creases the counts by about 0.5-1 dex, depending on whether 
or not Hq and baryonic acoustics oscillation constraints are 
included as well (for which the differences from the WMAP 
alone results are smaller). 

Most of the features arise in haloes with masses between 
10^ — 10® AIq . More massive haloes produce fewer absorption 
features both because they are fewer in number and because 
of their higher gas temperatures, which act to reduce the op- 
tical depth. The full mass range of the haloes that dominates 
the counts is beyond the reach of current cosmological nu- 
merical simulations. The results presented here indicate the 
mass resolution required to capture the statistics at a given 
equivalent width limit. Because higher mass haloes do not 
contribute much to systems with observed equivalent widths 
between 0.1 — 0.3 kHz, the statistics are fairly insensitive to 
the role of star formation in limiting the upper mass range 
of the haloes: results with and without molecular hydrogen 
formation yield similar results over this equivalent width 
range. 

The observed FWHM widths of the features are typi- 



© 2011 RAS, MNRAS 417, 1480-1509 



The micro -structure of the intergalactic medium I: the 21cm signature from dynamical minihaloes 27 



cally 4 kHz, so that the absorption features would be readily 
resolved in frequency by currently existing or planned radio 
facilities capable of detecting a cosmological 21cm signal, 
given a sufficiently bright background source to build up the 
required signal to noise ratio. Illustrative cases for detections 
by SKA and LOFAR are presented, for both individual ab- 
sorption lines and the ensemble spectral fluctuations over 
broad frequency ranges they will induce. It is shown that 
an observation able to detect the ensemble fluctuations will 
also detect the stronger individual absorption features, given 
adequate spectral resolution. 

The haloes have typical angular diameters of 0.05- 
0.25 arcsec ai z — 10. Whilst the radio core of a quasar 
is smaller than this, it is unclear that the dominant radio- 
emitting region of a background bright radio galaxy would 
be. The lack of full coverage of the emitting region would 
reduce the number of absorbers detected. It could alter the 
equivalent width distribution as well. If the larger haloes giv- 
ing rise to the larger equivalent width absorbers completely 
covered the emitting region while the smaller did not, then 
there would be a large relative reduction in the expected 
number of weaker equivalent width absorbers. The equiva- 
lent width distribution could thus serve as a probe of the 
size of the emitting regions of radio galaxies. 

In a scenario allowing for the production of a metagalac- 
tic UV radiation field, a Lya photon scattering rate match- 
ing the thermalization rate, required to decouple the spin 
temperature from the CMB temperature and begin coupling 
it to that of the IGM through the Wouthuysen-Field effect, 
produces absorption features with broad absorption wings. 
The equivalent width is computed relative to the wings at 
an observed frequency offset of 4 kHz. In some cases, mock 
emission lines result relative to the broad absorption wings. 

In the presence of Lya photon scattering at a rate 
matching or exceeding the thermalization rate Pth, the num- 
ber of absorption features is greatly enhanced over the case 
with no scattering, by more than two orders of magnitude 
at z = 10 for systems with w°^^ > 0.15 kHz. Increas- 
ing the Lya photon scattering rate to lOPth, so that the 
spin temperature is strongly coupled to the IGM temper- 
ature everywhere, increases the number of systems with 
10°^= > 0.25 kHz by three orders of magnitude over the case 
with no scattering. 

Adding a moderate amount of heat to the gas, how- 
ever, substantially suppresses the number of systems. A sud- 
den temperature boost by ten times the CMB temperature 
at 2 = 10 reduces the number of absorption systems with 
vjZ^" > 0.2 kHz by about 1 dex in the absence of Lya 
photon scattering from the case with no temperature boost. 
Allowing for a Lya photon scattering rate ten times the 
thermalization rate reduces the number of systems by close 
to a further order of magnitude. For gradual heating up to 
ten times the CMB temperature, adiabatic compressional 
heating dominates over shock heating within the core, pre- 
venting the formation of the high gas density central peak 
produced in the case without external heating. As a result, 
absorption features with wf^^ > 0.13 kHz are eliminated 
altogether. 

Heating of the IGM is thus degenerate with the sup- 
pression of small scale power in terms of the number of de- 
tectable absorption features against a bright background ra- 
dio source. On the other hand, probing an unhealed IGM, 



even if only the remaining cold patches as heating sources 
begin to turn on, would result in an unmistakable signal of 
a cold IGM because of the very large number of absorption 
systems that result. 

For a scenario with no large-scale galactic feedback, the 
minihaloes will emit relative to the CMB. The individual fea- 
tures have typical observed FWHM widths of 4 kHz, but are 
unlikely to be resolvable at the signal-to-noise levels achiev- 
able in the foreseable future. Their collective signal, how- 
ever, may be detectable. Because the strength of the signal 
of an individual halo increases with the mass of the halo, 
the overall signal is very sensitive to the assumed maximum 
halo mass. For models including molecular hydrogen forma- 
tion, the brightness temperature differential at z — 15 is 
STb = Tb — Tcmb — 0.2 mK, increasing to 0.6 mK at z = 8. 
If molecular hydrogen formation is suppressed, the increase 
in the upper halo mass raises the signal to 0.4 mK at z = 15 
and 1 mK at 2: = 8. At these levels, the minihalo signal is 
comparable to that of the diffuse IGM. 

Allowing for a running spectral index constrained 
jointly by the WMAP and ACT data suppresses the bright- 
ness temperature by about a factor of 30 at z = 15, and a 
factor of 5 at 2: = 8, without adding the constraints from 
_ffo and baryonic acoustic oscillation measurements. Adding 
these constraints reduces the suppression factors to 8 at 
z = 15 and 2-3 at z = 8. 

The 21cm signature against the CMB is extremely sen- 
sitive to large-scale galactic feedback. Allowing for a meta- 
galactic UV background that produces a Lya scattering 
rate as small as one percent of the thermalization rate re- 
sults in an overall absorption signal from the gas within a 
few times the turn-around radius of the haloes. The overall 
signature from the minihaloes, including their environment, 
would be one of absorption. 

The absorption signal would weaken as the IGM were 
heated. Once the IGM temperature exceeds the CMB tem- 
perature, however, the opposite effect occurs. Even with- 
out Lya scattering, the emission signal from the gas be- 
yond the core of the halo swamps the signal from within the 
minihalo core. Allowing for Lya scattering further increases 
the strength of emission from the surrounding gas compared 
with that of the core. 

Thus, a minimal amount of galactic feedback, either 
through the production of Lya photons or through heating, 
results in a signal against the CMB from the diffuse IGM 
that overwhelms that of the minihaloes. The direct mea- 
surement of a minihalo signal prior to any galactic feedback 
would require an absolute calibration of the brightness tem- 
perature of a radio telescope to an accuracy of better than 
STb /TcMB — 10~* for a detection at 2 = 15. More practi- 
cal would be an experiment designed to measure the differ- 
ential brightness temperature between the minihaloes and 
the CMB by differencing the signals from neutral and ion- 
ized patches, either in frequency or in angle on the sky. The 
circumstances for such a detection, however, would appear 
to be highly unlikely, perhaps achieved only in the earliest 
stages of reionisation, since the reionisation of a large frac- 
tion of the IGM would almost certainly be accompanied by 
galaxy or AGN produced Lya photons, and possibly heat- 
ing as well, that w ould be prevalent throughout the IGM 
iMadau et al.lll997l ). 
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APPENDIX A: SPHERICAL TOPHAT HALO 
TEMPERATURE 



The temperature (in K) from RECFAST l|Seager et alfcOOOl 'l 
is fit to better than f 0% accuracy by 

1.95 



Tk 



0.023(1 + z) 
0.146(1 + 2)^-^° 
1.46(1 + 2)'-''' 

TCMB 



; 6 < 2 < 60 
; 60 < 2 < 300 
300 < 2 < 500 
z > 500. (Al) 



The residual ionisation fraction We/wH.Tot is fit to better 
than 5% accuracy by 

X 10^ = 1.8 + 0.00176(1 +2)'-^ ;6<2<30 



?^H,Tot 



2.0 + 0.00082(1 +2)'-"" 
2.0 + 0.00366(1 +2)'-' 
0.0124(1 + 2) 



; 30 < 2 < 100 
; 100 < 2 < 300 
; 300 < 2. (A2) 



The spherical collapse model forms a two parame- 
ter family of solutions for a growing spherical perturba- 
tion. A convenient choice for the parameters is the to- 
tal halo mass and epoch of collapse. For the tophat col- 
lapse of a halo of mass AI and comoving radius ro, with 
M = (4/3)7rpAf (2)[ro/(l+2)]-^ where pm{z) is the total mass 
density at high redshifts z when vacuum energy contributes 
negligibly, the infall velocity win is given by 



1 2 GM 



ZtxV'^ GM 



(l + 2c), 



(A3) 



where Zc is the collapse redshift (|Meiksin|[2009h . The cor- 
responding post-shock gas temperature in the strong shock 
limit, taking the shock velocity to be iiin, is given by (for a 
monatomic gas) T^h = (3/4)^rviriai, where (3/2)fcBrviriai = 
GMffi/r-v, or 



1/3 



(1 + Zc) 



9 f2Ti'^ 
~3~ 



1330(1 + 2.) MiJ^K 



GMfh 



72.1(1 + 2c 



ro/kpc 

M 
106 Mo 



2/3 



(A4) 



where m is the mean mass per particle, assumed to be neu- 
tral primordial hydrogen and helium. This is slightly below 
the nominal mean "binding temperature" of the gas, defined 
by (3/2)fcBTbind/m = E^i,,JMb ^ i3/5)GM/ry for a bary- 
onic mass in the halo of Mb, and approximating the halo as 
having uniform density. The temperature may be expressed 
as 

_ 2GMm_16 

-'bind - c- , - -.r-fsh 



76.9(1 + 2c) 



M 



106 Mp 



2/3 



K. 



(A5) 



These temperatures may be compared with the temperature 
resulting from the adiabatic compression of the gas in the 
collapsed halo. For an increase in the gas density by the fac- 



tor Ph/PM = fc 



ISiv , the temperature will increase by 



the factor T^a/Tk = {ph/pK 



N2/3 



= f~^. Using Eq. dXTJ, the 



resulting adiabatic compression temperature for 6 < 2c < 60 



Tad 



: 0.73(1 + 2c)'-®'^K. 



(A6) 



When this temperature exceeds Tsh, the gas in the collapsed 
halo will contract adiabatically rather than shock. This will 
occur for halo masses smaller than a minimal halo shock 



Me, 



1020(1 + 2c)' 



'Mr; 



(A7) 



If the IGM is heated, Msh will increase. For example, scaling 
the unperturbed IGM temperature by the CMB tempera- 
ture according to /Tcmb(2) results in the minimum shock 
mass Msh ^ 1.3 X ICpf^ Mq, independent of redshift. 



APPENDIX B: 
FUNCTION 



THE MINIHALO MASS 



The density distribution of minihaloes is uncertain. The res- 
olution requirements for an A'^-body simulation determina- 
tion are currently too prohibitive to cover the full range of 
required scales, although progress has been made from very 
large-scale simulations in estimating the halo number den- 
sity. A minimum resolution requirement may be estimated 
as follows. For a simulation to represent a fair sample of the 
universe, the rms density fluctuation across it should be very 
small. At 2 = 10, o" < 0.1 requires a box with a comoving 
side exceeding 15/i~' Mpc. From Eq. (I26|l . the Jeans mass 
at 2 = 10 is Mj ~ 7000Mq. For h = 0.70, these values cor- 
respond to a minimum of 5 X 10'° Jeans mass particles, and 
an even greater amount to resolve Jeans mass haloes. 

Nonetheless, substantial progress has been made to- 
wards quantifying the abundances of increasingly low mass 
haloes. Several such fitting formulas have been presented in 
the literature, but not down to th e halo masses required. The 
simulations of lReed et al.l (|2007l ). however, have reached to 
M > 3 X IQt'Mq at 2 > 10. They find the halo abundances 
are described to 4% rms accuracy by the fitting formula 



dn _ Q.mpcvit{Q) 
dM ^ 7ri/2M2 



exp(-t^)_H'(nnoff,t) 



dt 



dlogM' 



(Bl) 



where 



H{nncH,t) = 0.542 [l + 0.901t" 



X exp 



0.236r 



''+0.6Gi(t) + 0.4G2(t) 
0.0369 



(B2) 



(nncff +3)2 

with Gi{t) = exp{-(logt - 0.576)V[2(0.6)2]}, G2(t) = 
exp{-(logt-0. 926)2/(2(0.2)2]} andn^s = 6dlog t/dlog M- 
3. 

Whilst it is straightforward to extrapolate Eq. (|B1|) to 
the lower masses required, it is unclear how reliable the re- 
sult would be, especially as the length scale moves more 
closely towards the flicker-noise limit for which rieft ~ —3. 
A useful guide may be offered by the density of peaks in 
a primo rdial Gaussian dark matter density field, as formu- 
lated bv iBardeen et all (fioSS ) (BBKS). The density of peaks 
filtered on a length scale 7?/ is given by 

(B3) 



dn 
dM 



f^m/5crit(0) 



where 7 — (y\l{(7Qai 



(^^) exp(-t2)G(7,2'/27t) 



dt 



dlogM' 



and r, — 3^^'^ai/a2 with 



a|(M,2) 



dlogkA^{k,z)W {k-Rf)k^ 



(B4) 



© 2011 RAS, MNRAS 417, 1480-1509 



The micro -structure of the intergalactic medium I: the 21cm signature from dynamical minihaloes 29 




10" 10' 
log,„ M (h-' M^) 




0.01 r 



0.001 



10« 10' 108 

og,„ M (h-' 



109 



Figure Bl. Dimensionless dark matter mass function 
{M^ / po)dn / dM , where po = f2mPcrit{0), at z = 20 (top 
panel) and z = 10 (bottom panel). Shown are the fit from 
iReed et al.l ||2007|) (soUd fines) and the BBKS approximation 
(dashed lines). 



W{k;Rf) is a filter function, and G is a function provided 
by BBKS. 

The uncertainty in relating the peaks in the primor- 
dial density field to collapsed haloes that form later cre- 
ates an uncertainty in the relation between the filtering 
scale R[ and halo mass M, with, for a Gaussian den- 
sity profile for examp le, M = (27r)'^/'^f2mpcrit(0)r5 , where 
rs = fsRf (Bond fc Mversll 19961) . The BBKS approximation 
requires fixing /a, the approach adopted here for its simplic- 
ity. More sophisticated approaches allow for the tidal forces 
that modify the evol ution of the primordial perturbations 
l|Bond fc Mverj|l996t ). which can be particularly severe for 
small haloes in the vicinity of larger scale overdensities. 

The value for fs is found by matching the BBKS ap- 
proximation to Eq. (|BH) over the mass range 3 x 10^ < 
M < 10® h~^MQ at z = 10 and 20. Choosing a Gaussian 
filter function W was found to require a slow variation of 
/a over the mass range. A more effective filter that better 
describes the structure of haloes is a tophat convolved with 
an exponential, W = WthVKcxp, where 

W-Tnik; Rf) = I [sm{kRf) - kRf cos{kRf)] , (B5) 



ikRfr 



and 



cxp (^) -^/) 



l + {OAkRff] (B6) 

(N. Dalai, Y. Lithwick, M. White, personal communication). 
Adopting 5c = 1.683 was still found to require a variable fs, 
but a constant value is possible choosing instead Sc — 1.8 
at z — 10, with fs = 1.1, and 5c = 1.5 at z — 20, with 
fa = 0.8. The resulting fits are shown in Figure [Bl] Better 
than 15% agreement is found between the BBKS fit and 
Eq. (fBlj) over the mass range 3 x 10^ < M < 10® Mq at 
z — 10, and better than 10% agreement over the mass range 



3 X 10^ < M < 10* H-^Mq at z = 20. The extrapolation of 
Eq. ([BlJ to Afj < Af < 3 X 10^ H'^Mq matches the BBKS 
predictions to better than ~ 5% at z = 10 and to better 
than ~ 15% at z = 20. This suggests either relation may 
be used to estimate the number of minihaloes. In this paper 
predictions are based on Eq. (|B1[) . 



APPENDIX C: MOLECULAR HYDROGEN 
FORMATION 

CI Reaction networks 

Gas-phase molecular hydrogen production occurs princi- 
pally through two catalytic processes, via the formation of 
the intermediaries H~ or Ht: 



H« + e- 



and 



H-+7, 
> H2 -I- e" 



H" + H+ ^ H+ + 7, 
H" -f- -5> H2 -f 



(CI) 
(C2) 

(C3) 
(C4) 



Because of the low binding energy of the extra elec- 
tron in H~ , the reverse of reaction Eq. (|C1|I induced by the 
Cosmic Microwave Background also becomes important at 
high redshifts. At high temperatures, the molecular hydro- 
gen may be destroyed through: 



H" + H2 ^ 3H°. 



(C5) 



Molecular hydrogen may also be destroyed through collisions 
with protons (H+ -I- H2 -I- H"), but this will generally 

be negligible compared with Eq. (|C5|I in a largely neutral 
gas. 

For primordial abundances, formation via H~ gener- 
ally dominates ( Palla et al. 1983; Lepp & Shull 1984|). Much 
more extensive lists of reactions relevant to H2 formation 
have been explore d, including reactions wit h deuterium, he- 
lium and metals (|Glover fc JappsenllioOTl ). but these gen- 
erally contribute negligibly to the overall abundance of H2 
created in a largely neutral diffuse medium, although they 
may play roles in the cooling and collapse of a molecular hy- 
drogen cloud. Many of the rates are still poorly determined, 
and this may affect the rate of cooling and collapse once 
initiated by a sufficient abundance of molecular hydrogen 
del over fc Abel 2008). Since only the formation of sufficient 
molecular hydrogen to drive cooling is considered here, only 
the H~ process is included in the halo collapse models. For 
the sake of completion, however, the reaction rates involving 

are listed as well. 



C2 Reaction rates 

The reaction chain of Eqs. (IC1|) and (|C2[I corresponds to the 
system 



dt 
dt 



= fcinjjorie — (fc2nH0 -I- fc5i)njj- , 
= fc2Wjjonj]- — fegJlHaWHO, 



(C6) 
(C7) 
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where rijjo, rijj- , riHa and Ue are the number densities of H", 
H~ , H2 and e~ , respectively. The reaction chain of Eqs. (|C3|) 
and (|C4|I corresponds to the system, including formation by 
H-, 



dt 

driHQ 
dt 



(C8) 



(C9) 



are the number densities of H"*" and 



where njj+ and n„+ 

H^, respectively. 

The reaction rate coefficients use d are summarised 
in Ta ble ICll following the labelling in iGlover fc JappsenI 
|2003). A few comments are required. The rate s Rl.l a and 
Rl.lb are based on the computation of lde Jon3 (|l972l ). who 
finds 



fci ~ 8.92 X 10" 
where 



'T 



-3/2 



I3{T), 



(CIO) 



/3(r) 



dk 



+ 0.0555) [exp (i^) - exp (^^)] ' 

where T4 = T/IO* K and k is the wavenumber of the electron 
in the reverse reaction R51 (photo-detachment). Ide Jonel 
l| 19721 ) tabulates values for fei over 10 < T < 15000 K. 
The function I3{T) may be approximated as follows. Defining 
P{T) = y(a;)exp(0.875a;), where 



rp-l 



y{x) = 



dk- 



' (A;2 + 0.0555) [exp(15.8(fc2 + 0.0555)a:) 
Differentiating gives 

15^8^0.875. 



Both R3.1 and R3 .2 are fits to the values tabulated in 
iRamaker fc Peekl l| 19761 ). The fit R3.1 matches the tabulated 
values to better than 3 per cent for 100 < T < 32000 K, 
decreasing in accuracy at larger values until 17 per cent too 
small at r = 5 X lO'^ K. Between 2 < T < 30 K the fit 
is as much as 50 per cent too large. The fit R3.2 matches 
the tabulated values to within 40 per cent over the range 
10 < r < 32000 K, and to better than 15 per cent over the 
range 500 < T < 8000 K. 

The rate R4.1 agrees well with the value kj ~ 6.4±1.2x 
10"^" cm"^s~^ measured bv lKarpas et ah I (1 19791), but exceeds 
the al ternative rate R4.2 determined bv lPrasad fc Huntress! 
(|l98(]|) by a factor of several. 

Rate R7 has had a wide range of estimates, as recog- 
nised bv ISavin et al] (|2004bl ). who argue for a similar but 
more definitive estimate than provided by iGalli fc Pallal 
(|l998l ), which is smaller by nearly an order of magnitude 
compared with the earlier cross-section rough estimate of 
Ide Jond (|1972| ). 

Rate R9.1 is a fit provided by IShapiro fc Kan3 (Il987l ) 
to the cross section values provided by Chernoff, HoUen- 
(Cll)bach fc McKee in a personal communication. It is for hy- 
drogen de nsities nn <C 10 cm" when T < lO"* K. The rate 
R9.2 from iDove fc Mandvl (|l986h accounts for quasi-bound 
states, and is valid for nn < 100 cm~^. The rate R9.3 applies 
in the low density limi t nu < 1 cm~^ (|Lepp fc Shulllfl983l : 
iMac Low fc Shulllligiel ). 

The computations for this paper generally used rates 
Rl.la, R2.1, R9.2 and R51, although comparison runs were 
made using some of the alternative rates. The results were 
not found very sensitive to the choices. 



.(C12) 



dy 
dx 



r(5 / 2) 0.8753: 3 

(15.8x)5/2 2ei-753:'*'l^ ' 2- 



1), 



(C13) 



where $(q, |,1) = E^=o(" + is a Lerch tran- 

scendental function. The resulting infinite series represen- 
tation for /3(r) is resummed as a second order Fade ap- 
proximant, which is the form provided in Table | C 1 J The 
expression agrees with the tabulation of Ide Jonel |l972l ) 
to better than 1 per cent for T < 3000 K, and to better 
than 4 per cent for T < 15000 K. In the limit T > 1, 
the rate approaches the asymptotic value fci ~ 8.92 x 
10-"r(5/2)C(3/2)/[3(15.8)^/2] ~ 1.64 x 10-i^cm3s-\ 
where denotes the Riemann function. This is the value 
given by the Fade approximant for T ~ 10^ K. An alterna- 
tive r ate, R1.2a and R1.2b, has been pr ovided bylWishartI 
l|l979h . which agrees reasonably well with lde Jond l| 19721) for 
T < 15000 K, and extends the range to higher temperatures. 
At T = 10^ K, it corresponds to ki ~ 6.34 x 10"^^ cm^ s"\ 
slightly higher than the value given by Rl.la at T = 
15000 K. For this reason Rl.la is fixed at the value at 
T = 15000 K (which agrees well with R1.2b at this tem- 
perature). The rate for R51 is based on Rl.la assuming the 
incident radiation is black body with a temperature T. 

The l aboratory meas u remen t of rate R2 at T = 
300 K by iMartinez et~aLl Hooi) of ^ 2.0 ± 0.6 x 
10~^ cm^ s~^ agrees v ery we ll with the theoretical value from 
iDalgarno fc Brownd (Il967l ) (R2.1), but is a lso consistent 
with the somewhat lower value estimated bv iLaunav et al.l 
(119911) . 



C3 Numerical integration scheme 

Despite its simplicity, the system of Eqs. HC6|) and HC7|) is 
not straightforward to solve numerically because of the dis- 
crepant timescales between the formation of H" and H2, 
rendering the system stiff. The solution to the system is 

nH-(t) = nH-,cq(t)+ [nH-(0)-nH-,cq(0)] 



X exp 

riHa (t) = riHa (0) exp 
-I- exp 



dt'{k2it')niia{t') + ksi{t')) 



; (C14) 



dt' kg{t')n^o (t') 
dt' kg{t')niio {t') 









X 


/ dt' exp 






Jo 





dt kg{t )njjo(t ) 
[k2{t')n^o{t')n^-{t')\ , 



(C15) 



H-,oq(*) = ki{t)n^Q{t)ne{t)/[k2{t)n^a{t) + k^^i{t)] is 



where n 

the equilibrium number density of njj- , which is rapidly 
achieved at high redshifts on the timescale l/{k2n^o ~\- kz,i). 

For solving the equations in a homogeneous expanding 
universe, it is useful to recast them as 



dXfi- 
dr 

dXB_2 

dr 



- fio — diixu2 ; 
hiXji- — d22a;H2 



(C16) 
(C17) 
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Table CI. Gas phase reaction rate coefBcients ki (cm^s ^). 



No. 


Reaction 


Rate coeiHcicnt 


Temperature range 


Rof. 


Rl.la 


HO 


■f e" ^ H" +7 


1 08 X 10-l*r4 4+26z+8z^ 
T/i — T/in* K 7 — Ta In 87^ 

J. 4: J. / A.\J J. V J J. 4- / \J . O ( o 


10 < T < 15000 K 


1 


Rl lb 








-L ^ lOUUU -iV 




R1.2a 






r\nv(—^ 7 Rd^ -U i/ifT) 7fi9 -1- mfC) 1 '^9'? — ()'^')7d'in\^\ 


± \ UUUU -TV 


2 


R1.2b 






(Jex(— 16.420 + ?i)^(0.1998 + «)^(— 0.005447 + 4.0415 x 10~'^ui^))); 
= logio ^ 


T > 6000 K 




R2.1 


hO 


+ H- ^H2+e- 


1.3 X 10-9 (T/IO^K)""-^ 


100 < T < 32000 K 


3 


R2.2a 






1.5 X lO^'J 


T < 300 K 


4 


R2.2b 






4.0 X 10-9T-017 


T > 300 K 




R3.1 


HO 


+ H+ -> H+ + 7 


2.25x10-20 [l+^^;^7«3^^^ij,,,] 


1 < T < 10^ K 


5 








R3.2 






dex [-19.38 + to(-1.523 + id(1.118 - 0.1269to))]; 


10 < T < 3.2 X lO* K 


5 






H+ ^ H2 + H+ 


^ = logio T 






R4.1 


H + 


5.8 x 10-10 


T < 10* K 


1 


R4.2 






1 X 10-10 




6 


R7.1 


H+ 


+ Ha ^ H+ + HO 


3 X 10-10 exp(-21050/T) 


T < 10* K 


7 


R7.2 






10-^/(tu) exp(-21237.15/T); 

f{u}) = —3.3232183 + «)(3. 3735382 + «)(—l. 4491368+ 
«!(0.34172805 + «;(-0.047813720+ 

u)(0.0039731542 + «)(-0. 00018171411 + 3.5311932 x I0~^w)))))); 
w = log(T) 


100 < T < 3 X 10* K 


8 


R9.1 


hO 


+ Ha ^ 3H0 


1.12 X 10-10 exp(-70350/r) 




9 


R9.2 






8.04 X 10-iir|-0i2 exp(-51790/T)/(l + 0.2130r4)3-5i2; 

T4 = r/io*K 


10^ < T < 10^ K 


10 


R9.3 






6.67 X 10-i2ri/2/ exp(l + 63590/T) 




11 


R51 


H- 


+ 7 -> HO + e- 


1.04xlO«T,^/^|i^i^exp(-l/.) 
1.62 X lO^Tl^ exp{-l/z) 
T4 = T/IO^ K, z = r4/0.875 


10 < T < 10* K 
T > 10*K 


1 



REFE RENC ES: (lllde Jo ne (1972); (2) Wishart (1979); C3)'Dal garno fc Brownj lll967l): (4) Launay et al.l lll99ll): (5) Ramaker fc 
lll976l) ; (6) [Prasad fc Huntresa al98a): ('7),Galli fc Palla (il99i); (8) ISavin et al.N2004blla h: (9) IShapiro fc Kand jigSTl ): (10) 

iDove fc Mandvl lll986l) ; (lli lMac Low fc Shulll lll986l) 



Peek! 



where x^^- = ^H-Z'T-HToti = "Ha/jiHToti with "-HTot in- 
dicating the total density of hydrogen nuclei in all species, 
/lo = kineX^o/H{a), with Xj^o = nno/nH.Tot, rfii = 
{k2njiO+k5i)/H{a), /ai = k2niio / H (a) , daa = k<)n^o/H{a), 
H{a) is the Hubble constant at epoch a = 1/(1 + z), and 
T = log a. The rapidity with which R- reaches its equi- 
librium value suggests the following second-order accurate 
scheme 



+ a;. 



n \ 
'^H-,eqj 



exp 



n + 1 



+ 




(C18) 



/ai^H- exp 



(C19) 



for advancing the abundance fractions from time step t" to 
r"+i =r" + At. 



C4 Effect of molecular hydrogen production on 
minihaloes 

In this section, the effect of molecular hydrogen production 
on the collapse of minihaloes is discussed. The large densi- 
ties in the cores of minihaloes result in rapid cooling either 
through molecular hydrogen or coUisional excitation of hy- 
drogen. In the presence of molecular hydrogen cooling, the 
minimum halo mass for forming stars is found to range from 
0.9 to 1.6 X 10'' Mq for 30 > z > 8, declining to 0.4 x lO'' Mq 
at z = 50. The corresponding range of virial temperatures 
is 3700 - 1600 K. 

If molecular hydrogen is dissociated by the metagalac- 
tic UV radiat i on field fro m the first star s, as is argued by 
iHaiman et al. I (|l997 a) and 'Haim an et al.1 (|2000l ). the mini- 
mum halo mass for cooling sufficiently to form stars increases 
to 3 X 10^ - 8 X 10*' Mq for 50 > z > 8, corresponding to 
virial temperatures of 14000 — 4 600 K. This range ex tends 
the estimate of IQ- ^- ^ K g iven bv lHaiman et al.l (|2000l ). 

Hai man et al.l (|200C ) argued that only very low UV 
mctagalactic background levels are required to dissociate 
the molecular hydrogen in haloes with masses below 10^ — 
10* Mq . The argument was based on hydrostatic minihalo 
models for which it was presumed that the initial molec- 
ular hydrogen content in a halo was negligible. In fact, as 
shown in Figure [S] substantial molecular hydrogen forma- 
tion begins well before collapse, already for overdensities of 
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Mhaio/ Me r (kpc) 



Figure CI. Mass in molecular hydrogen within the turn-around 
radius, as a function of total halo mass, at collapse epochs Zc = 8 
(triangles), 15 (squares) and 30 (circles). Open symbols indicate 
haloes that formed more than 1000 Mq of stars by the time of 
their collapse. 



only a few. A column density of A^H2 — 5 x 10^* cm^^ is 
required for self-shieldi ng in a 1000 K hal o from an exter- 
nal UV radiation field jHaiman et al.ll200Gl ). This is compa- 
ra ble to, but ab o ut a factor two larger, than the estimate 
of lde Jong et al.l (|l980l ) for the required column density for 
unit optical depth due to Lyman band absorption alone, so 
that the optical depths shown below may be somewhat un- 
derestimated. In the model shown, for which the collapse 
epoch is Zc = 20, the molecular hydrogen column density 
reaches the shelf-shielding value hy z — 30. Allowing for 
internal motions, howe ver, may require highe r column den- 
sities for self-shielding (| Glover fc BrandllioOll ). 

The fate of early haloes in which molecular hydrogen 
was formed was also neglected in past computations. These 
haloes will merge into larger ones, pre-enriching them in 
molecular hydrogen. The amount of molecular hydrogen 
within the turn-around radius of non-pre-enriched haloes is 
shown at the time of their collapse in Figure [CT] This is the 
amount of molecular hydrogen formed in the diffuse gas re- 
tained in the halo, as distinct from molecular hydrogen that 
would form in the gas that becomes thermally unstable and 
is removed from the halo at the rate p, given by Eq. (|27p . 
presumed to form stars. The open symbols correspond to the 
molecular hydrogen mass within haloes able to form 10'^ Mq 
of stars by the time they collapse. Only a small amount of 
molecular hydrogen need be present to reduce the cooling 
time sufficiently for the stars to form. 

The earlier the collapse epoch, the larger the mass in 
molecular hydrogen for a given halo mass. Approximating 
the optical depth to Lyman- Werner radiation as tl-w — 
X 10"cm~2)(1000 K/r)^/^ Figure [C2l shows that 
the molecular hydrogen would be self-shielding. A merger 
event could further concentrate the gas towards the newly 



Figure C2. Optical depth of H2 to Lyman- Werner photons in 
haloes collapsing at Zc = 30 as a function of comoving radius. 
Shown for haloes with initial total masses log(M/M0) = 5.4 
(short-dashed line), 5.6 (long-dashed line) and 5.8 (solid line). 
The optical depth is computed inwards from the turn-around ra- 
dius of the halo. 

formed minihalo centre. Mergers of molecular hydrogen en- 
riched haloes could thus build a substantial reservoir of 
shielf-shielded molecular hydrogen clouds within the cores 
of the minihaloes by the time a substantial UV background 
develops, lowering the minimum halo mass required for star 
formation. The topic is worth further investigation with fully 
3D computations, including a more complete assessment of 
the degree to which external Lyman- Werner photons will 
penetrate the halo core in the presence of velocity gradients 
within the halo. 

As pre-enrichment by mergers was not considered in the 
models computed here, the minimum halo mass for forming 
a sufficient number of stars to disrupt the geis in the halo 
may be even smaller than the models suggest. This would 
further reduce the numbers of high equivalent width absorp- 
tion systems along the line of sight to a bright background 
radio source and weaken the emission signature from mini- 
haloes against the CMB. 
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